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Abstract. We prove that for evolution problems with normally hyperbolic trapping in 

phase space, correlations decay exponentially in time. Normal hyperbolic trapping means 
that the trapped set is smooth and symplectic and that the flow is hyperbolic in directions 
transversal to it. Flows with this structure include contact Anosov flows [21], [37], classi- 
cal flows in molecular dynamics [25], [26], and null geodesic flows for black holes metrics 
[17], [18], [42]. The decay of correlations is a consequence of the existence of resonance free 
strips for Green's functions (cut-off resolvents) and polynomial bounds on the growth of 
those functions in the semiclassical parameter. 



1. Statement of results 

1.1. Introduction. We prove the existence of resonance free strips for general semiclassical 
problems with normally hyperbolic trapped sets. The width of the strip is related to certain 
Lyapunov exponents and, for the spectral parameter in that strip, the Green's function (cut- 
off resolvent) is polynomially bounded. Such estimates are closely related to exponential 
decay of correlations in classical dynamics and in scattering problems. The framework to 
which our result applies covers both settings. 

To illustrate the results in a simple setting consider 

(1.1) P = -h^A + V{x), VeC^{W]R). 

The classical flow cpt '■ (2^(0),C(0)) h-?- {x{t),^{t)) is obtained by solving Newton's equations 
x'(t)(t) = 2^(t), ^'(t) = — Vl^(x(t)). The trapped at energy E, Ke, is defined as the set of 

(x,^) such that p{x, ^) =^ ^'^ + V{x) = E and ipt{x,^) y^- oo, t — )■ ±oo. 

The flow ipt is said to be normally hyperbolic near energy E, if for some S > 0, 

def 1^ j^^^ ^ smooth symplectic manifold, and 

(1.2) |E-E'|<5 

the flow ipt is hyperbolic in the directions transversal to , 

see (1.17) below for a precise definition, and [26] for physical motivation for considering 
such dynamical setting. A simplest consequence of Theorems 2 and 6 is the following result 
about decay of correlations. 
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Theorem 1. Suppose that P is given by (1.1) and that (1.2) holds, that is the classical 
flow is normally hyperbolic near energy E. Then for G C^{{E — 5/2, E -\-5/2)), and any 
f,ge L2(M"), with 11/11^2 = ll^ll^. = 1, supp/,supp5 C B{0,R), 

(1.3) |(e-^*^/V(P)/,^)L^(M")| < ^^^^^e-'^' + Cn,Nh^, t > 0, 

for any 7 < Aq and for all N. Here Xq and Cq are the same as in (1.18) and Cr, Cr^n (^re 
constants depending on R and on R and N, respectively. 

This means that the correlations decay rapidly in the semi classical limit: we start with 
a state localized in space (the support condition) and energy, ■ip{P)f, propagate it, and 
test it against another spatially localized state g. The estimate (1.3) is a consequence of 
the existence of a band without scattering resonances and estimates on cut-off resolvent 
given in Theorem 2. When there is no trapping, that is when Ke = 0, then the right hand 
side in (1.3) can be replaced by 0{{h/t)°°), provided that t > Te, for some Te - see for 
instance [30, Lemma 4.2]. On the other hand when strong trapping is present, for instance 
when the potential has an interaction region separated from infinity by a barrier, then the 
correlation does not decay - see [30] and references given there. 

More interesting quantitative results can be obtained for the wave equation or for decay 
of classical correlations: see §1.2 for motivation and [42, Theorem 3] and Theorem 5 below 
for examples. When the outgoing and incoming sets at energy E, 

r| = {{x, : P{x, = E, iptix, ^) oo,t too}, 

are sufficiently regular, Theorem 1 and Theorem 2 below (without the specific constant 
Ao) are already a consequence of earlier work by Wunsch-Zworski [42, Theorem 2], and, in 
the case of closed trajectories, Christianson [9], [10]. For a survey of other recent results on 
resolvent estimates in the presence of weak trapping we refer to [40] . 

When normal hyperbolicity is strengthened to r-normal hyperbolicity for large r (which 
implies that arc manifolds) and provided a certain pinching condition on Lyapunov 
exponents is satisfied, much stronger results have been obtained by Dyatlov [19]. In par- 
ticular, [19] provides an asymptotic counting law for scattering resonances below the band 
without resonances given in Theorems 2, 4 and 6. It shows the optimahty of the size of 
the band in a large range of settings, for instance, for perturbations of Kerr-de Sitter black 
holes. 

Similar results on asymptotic counting laws in strips have been proved by Faurc-Tsujii 
in the case of Anosov diffeomorphisms [22], and recently announced in the case of contact 
Anosov flows [23] . In the latter situation, described in Theorem 4 below, the trapped set is 
a normally hyperbolic smooth symplectic manifold, but the dependence of the stable and 
unstable subspaces on points on the trapped set is typically highly nonsmooth. 

In this paper we do not assume any regularity on F^ and provide a quantitative estimate 
on the resonance free strip. For operators with analytic coefficients this result was already 
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obtained by Gerard-Sjostrand [25] with even weaker assumptions on . A new compo- 
nent here, aside from dropping the analyticity assumption, is the bound on the Green's 
function/resolvent that aUows apphcations to the decay of correlations. 

The proof is given first for an operator with a complex absorbing potential. This al- 
lows very general assumptions which can then be specialized to scattering and dynamical 
applications. 

Finally we comment on the comparison between the resonance free regions in this paper 
and the results of [32] where an existence of a pole free strip was given for hyperbolic trapped 
sets, provided a certain pressure condition was satisfied. In the setting of [32] the trapped 
set is typically very irregular but, the assumptions of [32] also include the situation where 
is a smooth symplectic submanifold, and the flow is hyperbolic both transversely to 
and on each Ke- In that case the resonance gap obtained in [32] involves a topological 
pressure associated with the full (that is, longitudinal and transverse) unstable Jacobian, 
namely 



where the supremum is taken over all flow-invariant probability measures on and H{ij) 

is the Kolmogorov-Sinai entropy of the measure // with respect to the flow. The bound is 
nontrivial only if this pressure is negative. In the case of mixing Anosov flows discussed in 
§9 the transverse and longitudinal unstable Jacobians are equal to each other; the above 
pressure is then equal to the pressure 7'( — log J'^), equivalent with the pressure V{— log J") 
of the Anosov flow, which is known to vanish [6, Proposition 4.4], and hence gives only a 
trivial bound. For this situation, our spectral bound (Theorem 4) is thus sharper than the 
pressure bound. On the other hand, one can construct examples where the longitudinal and 
transverse unstable Jacobians are independent of one another, and such that the pressure 
(1.4) is more negative — hence sharper — than the value — Aq given in (1.19), which may 



Notation. We use the following notation g = Ok{f)v means that \\g\\v < Ckf where the 
norm (or any seminorm) is in the space V, an the Ck depends on k. When either k or V 
are dropped then the constant is universal or the estimate is scalar, respectively. When 
F — Ok{f)v-^w then the operator F :V has its norm bounded by Ckf. 

1.2. Motivation. To motivate the problem we consider the following elementary example. 

Let X = M and P = —dl- A wave evolution is given by U(t) =^ sm{-\fPt) / \fP . Then for 
/, 5^ e (M) and any time t e M we define the wave correlation function as 





(1.5) 
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In this 1-dimensional setting, the correlation function becomes very simple for large times. 
Indeed, for a certain T > depending on the support of / and g, it satisfies 

Vt>T, C{f,gm = l f fix) f g{x)dx 

This particular behaviour is due to the fact that the resolvent of P, 

R(X) (p - : L^(R) ^ L\R), ImA > 0, 

continues meromorphically to C in A as an operator L^^^^^ Lf^^^ and has a pole at A = 0. 
In this basic case we see this from an explicit formula, 

[R{X)f]{x)^^l^e^'\^-y\f{y)dy. 

More generally, we can consider P = -dl + V{x),V e L^{R), with V > 0, for simplicity 
With the same definition of U (t) we now have the Lax-Phillips expansion generalizing (1.5) 

(1.6) C{f,g){t)^ f U{t)fgdx^ J] e"^^^* / fujdx f guj dx + 0{e-'''), 

Jr imA,>-A Jr 

where Xj are the poles of the meromorphic continuation of -R(A) = (P — A^)~^ (for simplicity 
assumed to be simple), and uj are solutions to (P — X'j)uj = satisfying Uj{x) = a.sgnxe*''*''^' 
for I a; I 3> 1. Since Uj are not in their normalization is a bit subtle: they appear in the 
residues of P(A) at A^ . 

The expansion (1.6) makes sense since the number of poles of P(A) with ImA > —A is 
finite for any A. If we define C{f,g) to be for i < 0, (1.6) gives (provided is not a pole 
ofP(A)), 

(1.7) C{f,g){-X)= Y^^: + ^(t)' Cj'^^ -i I fujdx I gujdx. 



The Lorentzians 



Im Xj I 



-2Im 



|A-Ajf Xj-y 

peak at A = Re Xj and are more pronounced for Im Xj small. This stronger response in the 
spectrum of correlations is one of the reasons for calling Xj (or A|) scattering resonances. 

In more general situations, to have a finite expansion of type (1.6), modulo some ex- 
ponentially decaying error 0{e~^*), we need to know that the number of poles of P(A) is 
finite in a strip ImA > —7. Hence exponential decay of correlations is closely related to 

resonance free strips. 

This elementary example is related through our approach to recent results of Dolgopyat 
[16], Liverani [29], and Tsujii [37] on the decay of correlations in classical dynamics. 
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Sum of Lorentzians corresponding to resonances 




Resonances obtained using splinepot([0, 50,20,80,0], [-2 -10 12]) 




Figure 1. Tlie effect of resonances on on the Fourier transform of corre- 
lations as described in (1.7). Tlie resonances are computed using tlie code 
scatpot .m [3]. 



Let X be a compact contact manifold of (odd) dimension ra, and let jt be an Anosov flow 
on X preserving the contact structure — see §9 for details. The standard example is the 
geodesic flow on the cosphere bundle X = S*M, where (M, g) is a smooth negatively curved 
Riemannian manifold. Let U{t) : C~(X) ^ C°°{X) be defined by U{t)f = Ytf = f ° It 
and let dx be the measure on X induced by the contact structure and normalized so that 
vol(X) = 1. The results of [16], [29], [37] show that, for any test functions f,g E C°°{X), 
the correlation function satisfies the following asymptotical behavior for large times: 

(1.8) C(/,^?)(t) = / [Uit)f]ix)gix)dx= [ fdx [ gdx + Oie-^'), t ^ cx) , 

Jx Jx Jx 

and the exponent F is independent of f,g. In other words, the Anosov flow is exponentially 
mixing with respect to the invariant measure dx. 
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Prom the microlocal point of view of Faure-Sjostrand [21], this result is related to a 
resonance free strip for the generator of the flow 7(. The resonances in this setting are 
called Pollicott-Ruelle resonances. 

In this paper we consider general semiclassical operators modeled on P given in (1.1), 
for which the classical flow has a normally hyperbolic trapped set. Schrodinger operators 
for which (1.2) holds appear in molecular dynamics - see the recent review [26] for an 
introduction and references. In particular, [26, Chapter 5] discusses the resonances in some 
model cases and the relation between the size of the resonance free strip and Lyapounov 
exponents. As reviewed in §9, the setting can be extended such as to include the generator 
of the Anosov flow of (1.8), namely the operator P{h) on X such that U{t) = 7^ = 
ex^{-itP/h). 

1.3. Assumptions and the result. The general result. Theorem 2, is proved for operators 
modified using a complex absorbing potential (CAP). Results about such operators can then 
be used for different problems using resolvent gluing techniques of Datchev-Vasy [13] - see 
Theorems 3 and 4. The assumptions on the manifold X, operator P, and the complex 
absorbing potential may seem unduly general, they are justified by the broad range of 
applications. 

Let X be a smooth compact manifold with a density dx and let 

P = P{x,hD) e ^'"(X), 

be an unbounded self-adjoint semiclassical pseudodifTerential operator on L'^{X,dx) (see 
§3.1 and [43, §14.2] for background and notations), with principal symbol p{x,^) indepen- 
dent of h. Let 

W = W{x,hD) e'^^X), 0<k<m, W>0, 
be another operator, also self-adjoint and with /i-indepcndent principal symbol w{x,$^), 
which we call a (generalized) complex absorbing potential (CAP). We should stress that 
W plays a purely auxiliary role and can be chosen quite freely. 

If the principal symbols p{x, ^) e S'^{T*X) and w{x, ^) e S^{T*X), we assume that, for 
some fixed Cq > and for any phase space point {x, ^) e T*X, 

- iw{x, 01 > {in Co -Co, 1 + w(x, > (O'/Co, 
ex.Tp{tHp){x, is defined for all i e R. 

Here, for ^ e T*X we have denoted = 1 + ||^||^ for some smoothly varying metric on 
X, X I— >■ II • 11^, and by Hp the Hamilton vector field of p. The map cxp(tiJp) : T*X — > T*X 
is the corresponding flow at time t. This flow will often be denoted by ipt, the Hamiltonian 
p{x,^) being clear from the context. 

We will need an additional smoothness assumption on w: 
(1.10) \d''w{x,0\ < CMx,0'-\ < 7 < I, 

when w{x, ^) < 1. This can be easily arranged and is invariant under changes of variables. 
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P = P-iW e 



the CAP-modified P. The condition (1.9) means that the CAP-modified P is classically 
elliptic and that for any fixed z & C 



Ke is compact and consists of points in p ^{E) which never reach the damping region 
{p e T*X : w{p) > 0} in backward or forward propagation by the flow (pf 

We illustrate this setup with two simple examples: 

Example 1. Suppose that Pq = -^^A + V,Ve C^{m";R), suppV ^ B{0,Ro). Define 
the torus X = M"/(3i?oZ)", and W e C~(X; [0, oo)), satisfying 

W{x)=0, xeB{(},Ro), W{x) = l, X e X\B{0,2Ro), d'^W = OiW"^'^), 

(here we identified the balls in M" with subsets of the torus). The last condition can 
be arranged by by taking W{x) = x(kP ~ Ro)i>{x) where x(a;) = exp(— x~"^)a;!|_, and 
-0 e C°°(X, (0, oo)) is suitably chosen. The power of W on the right hand side can be any 
number greater than |. 

Because of the support properties of V, P =^ — /i^A -\-V & ^^(X) and P — iW satisfy 
all the properties above. The trapped set Ke can be identified with a subset of T^^^q 
and is then equal to the trapped set of scattering theory: 



Example 2. Suppose that X is a compact manifold with a volume form dx and a vector 
field S generating a volume preserving fiow {C-dx = 0). Then P = —ihS is a self-adjoint 
operator on L'^{X, dx), and the corresponding propagator exp{—itP/h) is the push-forward 
of the fiow 7t = exp(tS) generated by S on functions / e L'^{X,dx): e:>q){—itP/h)f — 

f o 1-t- 

To define the CAP in this setting we choose a Riemannian metric g on X, and a function 



If Ag is the corresponding Laplacian on X, we set W{x, hD) — f{—h'^Ag). 

Then the operator P — iW G \E'^(X) satisfies the assumptions above. The principal 
symbols p{x,^) — CC^x), "wix,^) — f{\\C\\x)j where the norm || • is associated with the 
metric g. 



{{x,0 : p{x,0-z=p{x,0-twix,0-z^O}<^T*X. 
We define the trapped set at energy E as 

(1.12) KE = {p=ix,0 : pep-\E), 99m(p)c^/;-1(0)}. 



Ke = {{x, e T*W : + V{x) = E, x{t) 7^00, t^ ±00} . 



(1.13) 



/eC°°(K,[0,oo)), \f^'\s)\<Cf{sy-\ 0<7<|, 
f-\0) = [-00, M] for some M > 0, f(s) ^ y/s, s> 2M. 
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The trapped set is given by the points which never enter the damping region: 

KE = {{x,OeT*X : a^,) = E, ||(7-t)4L <M, Vt G M}. 

At tills stage the trapped set seems to depend on the choice of M. Below we will be 
concerned with exp(tS) being an Anosov flow, in which case this explicit dependence will 
disappear. 

Returning to general considerations we also deflne 

(1.14) ""^^ U Ke, 

\E\<5 

which is a compact subset to T*X and assume that 

(1.15) dp\KS^ 0, n WFh{W) = 0. 

The first assumption implies that for \E\ < 6, the energy shell p~^{E) is a smooth hyper- 
surface close to w^^{0). The second assumption is motivated by applications: in scattering 
problems the CAP is introduced only away from the trapped set. This second assumption 
also means that the trapped set is unchanged upon small modifications of WF/i(VF). 

We now make the following normal hyperbolicity assumption on K^: 

(1.16) is a smooth symplectic submanifold of T*X, 
and there exists a continuous distribution of linear subspaces 

3 Ef cTp{T*X), 

invariant under the fiow, 

yt e M, icp,)^Ef = 
and satisfying, for some A > 0, C > and any point p G K^, 



T,K' n E^ = E; n E; = {0} , dim Ef = dt, T,{T*X) = TpK' © ® E; 



^^YJ^ f 'p ^p''^p L-J 5 — "--p ■-'t, -pv- ^P^- ^^p 

V' ^^>u, m^^tKPM\^^^{p)^ue -\\v\\p. 



yv G Ef, yt > 0, |M(^=Ft(p)^ll^^*(p) < Ce-^*' 

Here p i— )■ || • ||p is any smoothly varying norm on Tp{T*X), p G K^. The choice of norm 

may affect C but not A. 

We can now state our main result. 

Theorem 2. Suppose that X is a smooth compact manifold and that P and W satisfy the 
assumptions above. If the trapped set given by (1.12), (1.14) is normally hyperbolic, in 
the sense that (1-16) and (1-17) hold, then for any e > there exists ho, cq, Ci, such that 
forO <h < ho, 

\\(P - iW - z)-%2^L. < C,h~'+'^'^'/'log{l/h), 
^ ' for z e [-6 + 6,6 -e]-ih[0,Xo/2-e], 
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where 



(1.19) 




Here det is taken using any fixed volume form on E^, the value of Xo > being independent 
of the choice of volume forms. 

This theorem will be proved in §6 after preparation in §§4,5. The bound log{l/ h)/h on 
the real axis is optimal as shown in [4] and [15]. Using the methods of [13] the estimate 
(1.18) almost immediately applies to the setting of scattering theory. As an example we 
present an application to scattering on asymptotically hyperbolic manifolds, which will be 
proved in §8: 

Theorem 3. Suppose (Y, g) is a conformally compact n-manifold with even power metric: 
Y is compact, dY — {x — 0}, dx\QY^ g — {dx^ + h)/x'^ where h is a smooth 2-tensor 
on Y with only even powers of x appearing in its Taylor expansion at x = 0. If the trapped 
set for the geodesic flow on Y is normally hyperbolic, then the following resolvent estimate 
holds: 

\\x'^{- A, - (n - 1)74 - XY'x'^^h.^L^ < Co^, A > 1. 

The next application is a rephrasing of a recent theorem of Tsujii [37]; it will be proved 
in §9. We take the point of view of Faure-Sjostrand [21], see also [12]. 

Theorem 4. Suppose X is a compact manifold and '-ft '■ X ^ X a contact Anosov flow 
on X. Let E be the vector field generating ipt, and P = —ihS the corresponding semiclas- 
sical operator, self-adjoint on L'^{X,dx) for dx the volume form derived from the contact 
structure. 

Define the minimal asymptotic unstable expansion rate 

(1.20) Ao =^ liminf \ inf log det {d'-it\EM) , 

with Eu{x) C T^X the unstable subspace of the flow at x. 

For any t > there exists a Hilbert space, Htc (see (9.10)), 

C°^{X)cHtG{X)cV\X) 

such that {P — z)~^ : Etc — ^ Hta is meromorphic for Im^; > —th. 

Then for any small e, 5 > 0, there exist Hq, cq > and Ci > such that, taking any 
t > Ao/2 and any < h < ho, 

(1.21) ||(P-z)-i^,^^^,^ <Ci/i-^+^°i-^/Mog(l//i), ze[S,6-']-ih[0,Xo/2-e]. 

We remark that under a pinching condition on the Lyapunov exponents, the recent results 
announced by Faure-Tsujii [23] provide a much more precise description of the spectrum of 
P = —ihS on HtG- the Ruelle-PoUicott resonances are localized in horizontal strips below 
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the real axis, and the number of resonances in each strip satisfies a Weyl's law asymptotics. 
That is analogous to the result proved by Dyatlov [19], which was motivated by quasinormal 
modes for black holes. 

Theorems 3 and 4 have applications to the decay of correlations, respectively for the 
wave equation and for contact Anosov flows. As an example we state a refinement of the 
decay of correlation result (1.8) of Dolgopyat [16], Liverani [29] and Tsujii [37]: 

Theorem 5. Suppose that ^yt : X ^ X is a contact Anosov flow on a compact manifold X 
(see §P.i for the definitions) and that Aq is given by (1.20). 

Then there exist a sequence of complex numbers, /ij, 

> Im II j > Im /ij+i, 

and of distributions Uj^k,Vj^k £ ^'{X), < k < Kj, such that, for any e > 0, there exists 
J(e) such that for any f,g & C°°{X), 



f{x)-ilg{x) dx = 

X 

(1.22) J{e) K, 

/ fdx / gdx + Y.Y.^'^~''''^^M)vM + OUe-'^'''^~'^'% 

''^ ''^ j=l k=l 

for t > 0. Here dx is the measure on X induced by the contact form and normalized so 
that vol{X) — 1, and u{f), u e T>'{X), f e C°°{X) denotes the distributional pairing. 

For information about microlocal structure of the distributions Uj^k and Vj^k the reader 
should consult [21]. Here we only mention that (with the standard wave front set of [27]) 

where E*, = \J^^^El{x), and El{x) C T*X is the annihilator of RS^ + E,{x) C T^X, 
• — u,s. The spaces E,{x) appear in the Anosov decomposition of the tangent space (9.2) 

2. Outline of the proof of Theorem 2 
The proof proceeds via the analysis of the propagator for the operator 

where the function ^(x,^; h) belongs to a certain exotic class of symbols. Our G which is 
closely related to the escape function constructed in [31], it depends on an additional small 
parameter, h, which will be chosen independently of h. 

For a large toj any fixed F > and e > 0, we can construct G so that, for some constant 
Co, the following holds uniformly m. Q < h < Hq, Q < h < Hq: 

G{p) = 0{\og{l/h)), G{p)-G{v^M)>-Co. peT*X, 
^^■^^ G{p)-G{^.to{p))>2T, pep-\[-d,S]), d{p,K')>{h/h)K w{p)<e, 
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where is any given distance function in T*X. 

The proof of Theorem 2 is based on the following estimate. For some eo > 0, take an 
operator A e such that WFh{A) C p-^{{-5,5)) n w-^([0,eo)). We will prove the 

norm estimate 

(2.2) \\eM-itPG/h)A\\mx)-,mx) < C e-'^^°-'^/^ , M,\ogl<t< Mlogl, 

h h 

As a result, for Im^; > — (Aq — 2e)/2, 

(2.3) {Pg -z)'- f e-^'^^^-'^'^Adt = (7 - g-^^^-^^"^)/'^) A = A - C»(e-^^)i2_i2 . 

Jo 

Hence, by taking T large enough and using the ellipticity of Pq — z away fromp~^((— 5, S))n 
w^^([0,eo)), we obtain 

{Pg - z)-' = Oih'')L2^Lh Im^ > -(Aq - 2e)/2 . 

Since e^'^'^ = 0{h~^)i,2^L2 from the growth condition on G, a polynomial bound for 
{P—iW—z)~^ follows. The more precise bound (1-18) follows from a semiclassical maximum 
principle. 

To prove the estimate (2.2) we proceed in a number of steps: 

Step 1. The most delicate part of the argument concerns the evolution near the trapped 
set. For some fixed R> 1, we introduce a cut-off function x G -Si supported in the set 



2 



{pep-\{-5,5)) : d{p,K')<2R{h/h)^. 

This cut-off is quantized into an operator == X^i^i hD). 
We then claim that 

for any e > and M > 0, there exists C > such that 

(2.4) 1 

uniformly for < i < M log — . 

h 

The proof of this bound is provided in §5. 

Step 2. For the weighted operator we obtain an improved estimate, now with a fixed large 
time to related to the construction of G, and for x which in addition satisfies 

X(p) = 1 for d{p,K') < R{h/h)K b(p)| < S/2. 

Using Egorov's theorem and (from (2.1)) the positivity of G — G o ip_-t;g on the set supp(l — 
x) n WF h{A), we get foUwoing the weighted estimate: 

(2.5) ||(l-x")e-^*°^«/M||<e-^ 
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When constructing the function G it is essential to choose V such that 

2 

We also show that 

(2.6) \\e-''°^°/^A\\ < e^^\ 

for a constant Cq independent of /i, h. Formally, these results follow from Egorov's theorem 
but care is needed as G is a symbol in an exotic class. To obtain (2.5) and (2.6) we proceed 
as in the proof of [31, Proposition 3.11]. This is done in §6. 

Step 3. The last step combines the two previous estimates, by decomposing 

In order to apply (2.4) we use the fact that 

where the symbols XG,i have the properties required in (2.4). A clever expansion of 
^^mtoPc/h^n jj^^Q terms involving Ug,± and an application of Steps 1 and 2 lead to the 
estimate (2.2) for t = nto. The argument is presented in §7. 

3. Preliminaries 

In this section we will briefly recall basic concepts of semiclassical quantization on man- 
ifolds with detailed references to previous papers. 

3.1. Semiclassical quantization. The semiclassical pseudodifferential operators on a 
compact manifold X are quantizations of functions belonging to the symbol classes S"^ 
modeled on symbol classes for R": 

5-(rM") = {ae C~(T*M" X (0,1],) : Va, /3 G N", \d:d^^aix,^;h)\ < Cafi{l + l^D^^-H^l} , 
see [43, §14.2.3]. The Wcyl quantization, which we informally write as 

S"XT*X) 3 a(,T,0 ^ a'"{x,hD) G ^"(X), 
maps symbols to pseudodifferential operators. It is modeled on the quantization on M": 
[a'^ujix) ^a'"{x,hD)u{x) = [Op^{a)u]{x) 




{27rhY 
The symbol map 

a : ^""{X) S"'{T*X)/hS"'-\T*X), 
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is well defined as an equivalence class and its kernel is h^'^~^{X) - see [43, Theorem 14.3]. 
If cr{A) has a representative independent of h we call that invariantly defined element of 

S"^{T*X) the principal symbol of A. 

Following [11] we define the class of compactly microlocalized operators 
^comp(^) def ^^^^^^ hD):ae {S^ n C^){T*X)} + 
These operators have well defined semiclassical wave front sets: 

see [11, §3.1] and [43, §8.4]. 

Let u — u{h), ||xi(/i)||L2 — 0{h~^) (for some fixed N) be a wavefunction microlocalized 
in a compact set in T*X. in the sense that for some A e ^^°™p, one has u — Au-\-Oc'^{h°°). 
The semiclassical wavefront set of u is then defined as: 

(3.2) WF;,(ii) = C{per*X : 3ae5°(r*X), a{x,0^l, \\a'"u\\L2^0{h°^)] . 
When A e *^°"^p(X) we also define 

WFh{I - A) := y WFh{B{I - A)), 

and note that WFh{B) n WF^(A) is defined for any B e *"*(X) as 

WFh(B) n WFh(A) := WF^(CS) n WFh(A), C E ^"""'p, WF^(7 - C) n WFh(A) = 0. 

Semiclassical Sobolev spaces, Hf^{X) are defined using the norms 

(3.3) \\u\\Hs^x) = \\iI-h'A,r/\h.^x), 

for some choice of Riemannian metric g on X (notice that H^{X) represents the same 
vector space as the usual Sobolev space i7*'(X)). 

3.2. Si calculus with two parameter. Another standard space of symbols ^'^(M^"), 

< S < 1/2, is defined by demanding that d"a = 0{h~^"'^^). The quantization procedure 
a (-> Op'^a gives well defined operators and Op^a o Op^b = Op^c with c G Ss- 

For 0<5<l/2we still have an pseudo differential calculus, with asymptotic expansions 
in powers of h. However, for 5 = 1/2 we are at the border of the uncertaintly principle, 
and there is no asymptotic calculus - see [43, §4.4.1]. To obtain an asymptotic calculus the 
standard Si spaces is replaced by a symbol space where a second asymptotic parameter is 
introduced: 

Then the quantization a i-> a"'(x, hD) e (R") is unitarily equivalent to 

(3.4) h^ay{x,hD) = Opl{a), a{p) = a{{h/hY^ p). a e ^(M^"), 
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— see [43, §§4.1.1,4.7.2]. Hence, we now have expansions in powers of h, as in the stan- 
dard calculus, with better properties (powers of {hh)^) when operators in and ^ are 
composed — see [35, Lemma 3.6]. 

For the case of manifolds we refer to [11, §5.1] which generalizes and clarifies the presen- 
tations in [35, §3.3] and [42, §3.2]. The basic space of symbols, and the only one needed 
here, is 

S'T^'iTX) ={a e Cr(T*X) -.V^-'-Vka^ 0{{h/h)--^), V k, 

2 

Vj e C^{T*X,T(T*X))} + h^S-^(T*X). 

The quantization procedure 

ST"'''{T*X) 3a^ Op^a) e ^'^^(X) 

2 2 

defines the class of operators ^^"^"^{X) modulo h°°^~°°{X), and the symbol map: 

2 

(3.5) a : ^T^iX) ST'''iT*X)/h^UsT'''{T*X). 

2 2 2 

The properties of the resulting calculus are listed in [11, Lemma 5.1] and we will refer to 
those results later on. 

When /i = 1 we use the notation 5*1°™^ (T*X) for symbols and denote by \I''i°™^(X) the 

2 2 

corresponding class of pseudodifferential operators. The symbol map 

a : *r"^P(X) ST''{TX)/h^ST''{T*X), 

2 2 2 

is still well defined but the operators in this class do not enjoy a proper symbol calculus in 
the sense that a{AB) cannot be related to a{A)a{B). However, when A e ^*f'™^(X) and 

2 

B e *(X) then a{AB) = a{A)a{B) +0{h^)s,{T''X) — see [35, Lemma 3.6] or [11, Lemma 
5.1]. 

3.3. Fourier integral operators. In this paper we will consider Fourier integral operators 
associated to canonical transformations. It will also be sufficient to consider operators which 
are compactly microlocalized as we will always work near p~^{[—2S, 25]) nty"^(0) which by 
assumption (1.9) is a compact subset oiT*X. 

Suppose that Yi,Y2 are two compact smooth manifolds {Yj — X or Yj — T"- in what 
follows) and that, Uj C T*Yj are open subsets. Let 

k:Ui^ C/2, F;, {{x, e, y, -77) : {x, = <y, v), iv, v) e U^} C ^ x T*Y,, 

be a symplectic transformation, for instance k = (pt, Ui = U2 = T*X. Here F^ is the graph 
of K, and ' denotes the twisting rj i->- —77. This follows the standard convention [27, Chapter 
25]. 
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Following [11, §5.2] we introduce the class of compactly microlocalized h-Fourier integral 
operator quantizing n, /^°™^(y2 x Yi, F'^). If T G /)j°™''(Y'2 x ^i; T'^) then it has the following 
properties: T = 0{l)L2(Yi)^L^Y2y, there exist Aj e *^°™p(Yj), WFh{Aj) ^ Uj such that 

A2T = T + 0{h°°)v'{Yi)-^C'^(Y2): TAi = T + C>(/i°°)x)'(Yi)-^C°°(Y2); 

for any Bj e 

TB,^CiT + hT,, a(Ci) ^ a(B^) o K-\ 
^ ' ^ B2T^ TC2 + hT2, a{C2) = a{B2) o Tj e ^"^^(^2 x Y,, T'J. 

The last statement is a form of Egorov theorem. 

When Bj e then an analogue of (3.6) still holds in a modified form 

2 

TBi = CiT + h^h'^DiTi, a(Ci) = a(Bi) o k-\ 

(3.7) B2T ^TC2 + h^h^T2D2, ^(Cs) = (7(^2) o 

2 

see Proposition 6.3 (applied with g = {)). 
An example is given by the operators 

(3.8) Ae-'^^'^, e-'*^'^A e r°"^P(X x X, T'^J, if A e *^°"^p(X). 

In §5 we will also need a local representation of elements of p°™'i' as oscillatory integrals 

- see [1],[20, §3.2] and references given there. If T G Jcomp(]R" x M'',!'^) is microlocalized 
to a sufficiently small neighbourhood k{U) xU G T*M" x T*M'^ ([43, 8.4.5]) then 

(3.9) Tu{x)^{2nh)-^ [ I ei'^^'''y''^^a{x,y,e)u{y)dyde + 0{K^)s\\u\\h-m , 

for any M. Here a G C~(M" x R''), ^{x,y,e) G C~(M2" x M'^), and near k{U) x U, the 
graph of K is given by 

T^ = {{{x,d^ilj{x,y,e)),{y,-dy'ilj{x,y,e)) : (x, 6*) G Cv-I' 

(3.10) c:/= {(x,|/,^) : de^(a;,|/,^) = 0,}, 
dx,y,0{dejip), j — 1, - ■ ■ ,k, are linearly independent. 

For given symplectic coordinates {x, ^) and {y, rj) in neighbourhoods of k{U) and ?7 respec- 
tively, such a representation exists with an extra variable of dimension k, where < k < n, 
and n + A; is equal to the rank of the projection 

assumed to be constant in the neighbourhood of k{U) x U - see for instance [43, Theorem 
2.14]. Since in (3.10) is the graph of a symplectomorphism it follows that for some 
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(3.11) 

For the use in §5 we record the following lemma, proved using standard arguments (see for 
instance [1]): 

Lemma 3.1. Suppose that T is given by (3.9) and that B e (M"). Then for any u E L'^ 
with \\u\\l2 — 1, 

BTu{x) = {27rhy^ [ [ eT^'^^'''y^'^^a{x,y,9)b{x,d:c^{x,y,9))u{y)dyd9 + 0{h'^U)L2, 

TBu{x) = {2TTh)-'^ I I eK^(^'^'^)a(a;, y, 9) b{y, -dyi){x, y, 9))u{y)dyd9 + 0{h^h^)L2, 
where b = (t{B). 

3.4. Fourier integral operators with operator valued symbols. In §5 we will also 

use a class of Fourier integral operators with operator valued symbols. We present what we 
need in an abstract form in this section. Only local aspects of the theory will be relevant 
to us and we opt for a direct presentation. 

Suppose that is a Hilbert space and Q is an (unbounded) self-adjoint operator with 

domain V <Z V.. We assume that Q : V ^ Ti is invertible and we put =^ Q'^H, for 
^ > 0. For £ < 0, we define T>^ as the completion of H with respect to the norm HQ^iiHw. 

We define the following class of operator valued symbols: 

(3.12) SsiM.'^'' xR\n,V), 
to consist of operator valued functions 

R" X R" X R'^ 9 (x, y, 9) ^ N(x, y,9) 
which satisfy the following estimates: 

(3.13) d:^^^,N{x, y, 9) = : ^ T^, 

for any multiindex a and £ e Z, uniformly in {x,y,9). We note that this class is closed 
under pointwise composition of the operators: if Nj e ^S^ then Nj defines a family of 
operators ^ V^, hence so does their product N1N2; the estimate (3.13) follows for the 
composition, since for \f3\ + I7I = \a\, 

d^Nid^N2 = 0{l)'j;)e+s\c,\^x>e+s(\c\-\p\)0{l)x>i+s(\a\-\i3\)^x>^ = 0{l)x>e+s\a\ ^x>'. 

Proposition 3.5 at the end of this section describes a class which will be used in §5. 

Suppose that ip satisfies (3.10) and (3.11). Wc can assume that -0 is defined on R^" x R'^. 
For N e Ss and a G C^(R^" x R'^) wc define the operator 

T : l2(r") ®n — > l2(r") ® n, 
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by 

(3.14) T{u^v) = {2TTh)-'^ I I eTi^^''^y^^^a{x,y,e){u{y)®N{x,y,e)v)dyde. 

This operator is well-defined since a is compactly supported, but to obtain a norm estimate 
which is uniform in h we need to assume that N & Sq: 

Lemma 3.2. Suppose that N e 5o(R^"+'=, 7^, X>) and that T is given by (3.14). Then 



(3.15) \\T\\L-2(Rn)^-H_,L2(^n)^ji = uiax — " /TF^^^ + 

where =^ {{x,y,6) : dgip = 0}, and is given by (3.11). 



IfNe cS^(M^"+^?^,^') then 

(3.16) T = 0(1) : L\W) (g) L\W) (g) V^, 
where rUn depends only on the dimension n. 

Proof. The estimate (3.15) follows from a standard argument based on considering T*T 
and from [43, Theorem 13.13]. The estimates (3.13) with 5 — Q and I — Q show that the 
operators can be treated just as scalar symbols. 

To obtain (3.16) we note that 

dly^9 {Q-'^Nix, y, 9)) = e»(l) :n^n, for \a\d < L. 

To obtain the norm estimate (3.15) we only need a finite number of derivatives, m„, de- 
pending only on the dimension. Taking L > rUnS, we can then apply (3.15) to the operator 
Q~^T, which gives the bound (3.16) for T. □ 

A special case of is given hj k — id. In that case we deal with pseudodifferential operators 
with operator valued symbols. The following lemma summarizes their basic properties: 

Lemma 3.3. Suppose that Nj e S5^{M?"), j = 1,2. Foru & andv e we define 

0^l{N,){u®v) j ei^-y^^ [A^,(^,OH u{y)dyd^. 

These operators extend to 

(3.17) 0^l{Nj) = 0{1) : L\W) ® V^^"^-^^ L\W) ® V^, 

and satisfy the following product formula: 
(3.18) 

0p;^(iVi)0p;^(iV2) = OpniViiV2)+^^, R = o(i) ■. l^(w')®v^+"'"^^'+^^^ L^{w')(g)V^. 

Here and in (3.17), £ is arbitray and rrin depends only on the dimension n. 
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Proof. When 5i — 82 — the proof is an immediate vector valued adaptation of the 
standard arguments presented in [43, §§4.4,4.5] where we note that only a finite number 
(depending on the dimension) of seminorms of symbols is needed. In general, (3.13) gives 

(3.19) d^^Q-^NjQ-^ = 0{l) : V\ \a\Sj <L + M, 

and the norm estimates (3.17) follows. To obtain the product formula we note that, using 
(3.19), it applies to Q^'^Ni and N2Q~^^ for M sufficiently large depending on n. Hence 

ovi{N,)Ovi{N2) = g^op]^(g-^iVi)OpniV2g-'")g^ 

= g'^0p]:;(g-^iViiV2g-^)g^ + g^c»(/i)L2^^.^i2^^.g^ 

which gives (3.18) provided m„(5i + ^2) > 2M. □ 

We can also factorize the operator T using the pseudodifferential operators described in 
Lemma 3.3, the proof being an adaptation of the standard argument. When 5" : L^(R") — >■ 
L2(K") we also write S for S ^ In : ^^(K'^) <^n^ L2(M") (g) n. 

Lemma 3.4. Suppose that T is given by (3.14) with N e .S^. Then 

T = rllOpHA^i) + hRi = Op;;'(A^2)Tll + hR2, 

where 
(3.20) 

Til e r°^P{W. X R",r;), rllti(x) = (27r/i)-^ / [ ei'^^'''y'^^a{x,y,e)u{y)dyde, 

N2{x,d^'i/j{x,y,e)) = Ni{y, -dytlj{x,y,e)) = N{x,y,e), {x.y.O) G 
Rj = C(l) : L\W) (g) V^"'"+^ L\W) ^ V^. 
Here, Nj e ^^(M" x R'',n,V), and 

Op'l^iNj) = 0{1) : L2(M") ® L^R"") ® V^. 

In our applications we will have 

(3.21) n^L^W^^dy), Q^-~eAy + f + l, 

so that T>^ are analogous to Sobolev spaces. In the rest of this section (as well in section 5), 
we will use the shorthand notations p|| = (x, y, 9) in order to shorten the expressions, and 
to differentiate between these variables and the "transversal variables" {y,fj). 

We consider a specific class of metaplectic operators: 

(3.22) Ar(p||)«(^) = (27r/.)-* / [ (deta?^g,„)^ei(^''i|(^'^)-<^'^'>^«(^')^^?/', 

jR'^t jRdt 

where {y, fj) is a real quadratic form in the variables y, rj, with coefficients depending 
on being in the class 5'(M^"'"'"'^), and the matrix of coefficients (?|,j^?p|| is assumed to be 
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uniformly non-degenerate for all p\\. The definition involves a choice of the branch of the 
square - see Remark 5.8 for further discussion of that. For any fixed these operators are 
unitary on "H (see for instance [43, Theorem 11.10]). 

The next proposition shows that this class fits nicely into our framework: 
Proposition 3.5. The operators N{p\\) given by (3.22) satisfy 

(3.23) d^N{p\\) = 0„,,(/i-l"l) : V\ 

for all I, That means that (3.13) holds with 5 — 1 (the loss in h is considered as dependence 
on a). 

IfX& ^(M^*) is fixed, A>1, and xa(«) = x(A"^«); then for any £ and k > 0, 

(3.24) xliV: hDy) = 0,(A2^) : ^ V'+\ 
so that 

(3.25) ) ( 

d^, (^(PIi)Xa(^, hDy)) = 0«,,(A^I"l/i-l"l) : V'. 

Proof. We see that dp^N{p\\) is an operator of the same form as (3.22) but with the ampli- 
tude multiplied by 

l^l<2|a| 

where < Hence to obtain (3.23), it is enough to prove that 

qV^7v(p||) ((^')''^(/^^^')''^Q-'-'"'^(y')) =0(lklk)H- 

Using the exact Egorov's theorem for metaplectic operators (see for instance [43, Theorem 
11.9]) we see that the left hand side is equal to 

N{p\\) {p'^piy'My') {K;Qy Q-'-\My')) , (dr^q,?,) ^ {y,dyq), 

q — qp^^ and where is a polynomial of degree less than or equal to \ j3\. Since < 2|q;|, 
the operator p'p{K*QYQ~^~\°'\ is bounded on % (see for instance [43, Theorem 8.10]) so 
the unitarity of N gives the boundedness in %. 

To obtain (3.24) we first note that xa ^ S'(M^*) uniformly in A > 1. Hence Q~^XaQ^ = 
(9(1)^^^, uniformly in A (again, see [43, Theorem 8.10]). This gives (3.24) for k = 0. For 
the general case we put Q\ = l+A~'^{{hDy)'^+y'^), and note that for any M, Q/^Xa — Xa,mj 
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where Xa,m G 5'(R^'^*) uniformly in A. Hence it is bounded on L^(R'^*) uniformly in A and 
h. We then write 

Q'^xliy, hDy) = Q^QtQixliv, hDy) 

= (1 + {hDyf + y^il + Pr^ih.Dyf + k-'e)-\lu{y. hDy) 
= + {hDyf + ff{K' + (/.Z),)^ + f)-\l,{y, hDy) 

completing the proof of (3.24). □ 

4. Classical dynamics 

In this section we will describe the consequences of the normal hyperbolicity assumption 

(1.16) , (1.17) needed in the proof of Theorem 2. 

4.1. Stable and unstable distributions. Let be the trapped set (1.14) and C 

TpX, p G K^, the distributions in (1.17). We recall our notation (pt =^ exptHp for the 
Hamiltonian flow generated by the function p{x,^). 

We start with a simple 
Lemma 4.1. If cu is the canonical symplectic form on T*X then 

(4.1) uj, 0, 

that is are isotropic. 

Without loss of generality we can assume that the distributions E^ satisfy 

(4.2) e;®e;^{t,kY: 

where V-^ denotes the symplectic orthogonal of V. 

Proof. The property (4.1) follows from the fact that ipt preserves the symplectic structure 
{ip;uj = uj). for X,YeE^, 

^p{^, y) = i^<p^t(p){{d^^t){p)x, d^^tip)y) ^ , t^+Qo 

To see that we can assume (4.2) we note that the distribution {{TpK^)-^, p e K^} is 
invariant by the flow: dipt{p) '■ TpK — )■ T^^(^p)K, and diptip) is a symplectic transformation. 

If Tip : Tp{T*X) — )• (TpK^)-^ is the symplectic projection, then Hdiptip) °dipt{p) = dipt{p) oTTp. 
This means that we may safely replace E^ with TTp{E^), without altering the properties 

(1.17) . □ 
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4.2. Construction of the escape function. To construct the escape function near the 
trapped set we need a lemma concerning invariant cones near . To define them we 
introduce a Riemannian metric on T*X and use the tubular neighbourhood theorem (see 
for instance [27, Appendix C.5]) to make the idetnifications 

neigh(ii'^) ~ N*K' n {(p,C) G T*{T*X) : ||C||p < 
(4.3) ~ {TK')^ n {(p, z) e T{T*X) : \\z\\p < 

~ {(m, z):meK\ ze R^''*, < ei} . 

Here {TK^)^ denotes the symplectic orthogonal of TK^ C Tks{T*X) C T{T*X). Since 
is symplectic, the symplectic form identifies (TK^)^ with the conormal bundle N*K^ . 
The norm || • ||p is a smoothly varying norm on Tp{T*X). We write dp{z^ z') = \\z — z'W^y and 
introduce a distance function d : neigh(ii''') x neigh(i^r'^) — )■ [0, oo) obtained by choosing a 
Riemannian metric on neig h{K^). We have ^((m,;^), (m, ;2')) ~ (im(2;, 2;') and the notation 
a b, here and below, means that there exists a constant C > 1 (independendent of other 
parameters) such that b/C < a < Cb. 

Assuming that are chosen so that (4.2) holds we can define (closed) cone fields by 
putting 

C7± = {z e {T,KY ■■ dp{z,Ei) < eMl, \\z\l < ej, 
C^"^ U C± C neigh(X^), 

where we used the identification (4.3). Since the maps p i-> E^ are continuous, are 
closed. 

The basic properties are given in the following 

Lemma 4.2. There exists to > and > such that, for every t > to there exists such 
that if one chooses ej < j — 1,2 in the definition of neigh{K^) and , then 



(4.5) p e C±, ip±t{p) e neigh(X^) =^ ip±t{p) e C±. 

In fact as stronger statement is true: for some constant Ai > and any t > to, 

(4.6) p, ip±t{p) e neigh(i^'^) d{ip±t{p), C^) < e~^''d{p, C±). 
Finally, 

(4.7) d{p,C^f + d{p,C-f^d{p,K'f. 

The conclusions (4.6) and (4.7) are similar to [31, Lemma 4.3] and [34, Lemma 5.2] but 
the proof does not use foliations by stable and unstable manifolds which seem different 
under our assumptions. 
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Proof. For p e neigh{K^) let (m, z),m e and z e R^''* ~ (T^X'^)-^ be local coordinates 
near p. Similarly let (m, z) be local coordinates near (pt{p) £ neigli(ir'') (by assumption in 

(4.5)). Then if for each m we put d±(fit{m) =^ d(fit{m)\(T^K^)-^, the map (ft can be written 
as, 

¥'t(m, ^) = (^t(m) + Oidkf ), rfxv^t(m)2; + Ot{\\z\\^)) 

^ {(pt{m,),d^ipt{rm)z + Ot{\\z\\'')), rm^m + Ot{\\zf). 

(Here we identify (TmK^)'^ with M^*^* and and consider d±Lpt{'m) ■ M^*^* M^"^', with similar 
identification near ift{p)- The norm || • || is now fixed in that neighbourhood.) 

Let z = Z-^- + be the decomposition of z corresponding to (Tmj^K^)-^ = © (we 
assumed without loss of generality that (4.2) holds). The continuity of p i— ?> E"^ and the 
definition of show that if ei is small enough depending on t (so that d{m, nii) = OtdlzW^) 
is small) 

(4.9) zeC+=^ Ik-IUi < 2e2|k+IUi. 

Since 

d±(pt{mi)z = ^ d±ipt{mi)z±, d±ipt{mi)z± e 
± 

and, by normal hyperbolicity, for some C > and Ai > one has for all positive times 

\\d±ipt{mi)z+\\ > ^6^^1*11^+11, 

\\d±^timi)z_\\ <Ce-^^'^\\z_\\. 

If 2; e C+, then this and (4.9) show 

\\d±Vt{mi)z^\\ < 2C^ e~*^'' e2\\d^(pt{mi)z+\\ . 

Let us take to such that 2C^ e-4Aito < 1/2. For t > to and ei small enough depending on t 
this shows that 

(4.11) zeC^ and ||2;|| < ei, \\d±ipt{mi)z\\ < ei ^ d±(pt{mi)z + Ctdl^lH G C+(^^), 

which in view of (4.8) proves (4.5) in the + case with the — case being essentially the same. 
To obtain (4.6) we note that for {m,z) G neigh(p, i^^), 

d{{m,z),C+) r.dm{z,C+) r.\\z4{l-lc+{z)), z^z+ + z^, z±eE^, 

where 1a is the characteristic function of a set A. (To see the first ~ wc need to show that 
CQdrn{z,C^) for some cq, which follows from an argument by condradiction 
using prc-compactness of .) 

We also observe that if diptip^i)^ G neigh(ii'^) then (4.11) gives, for ei small enough 
depending on t, 

1- V, id^^t{m,)z + Ot{\\z\\^))<l-lc^{z). 
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Hence, using (4.8) and (4.10), writing z — z^ + z^ as before, and taking ei sufficiently small 
depending on t >to, 

d{cpt{m,z),C+) ~ d^^{rm){d±^t{mi)z + Ot{\\zf),C+^^^^) 

~ \\d^^t{m,)z.\\ (1 + 0,{\\z.\\)) (1 - V, id±M^i)z + 0,{\\z\\f)) 

<Ce-''^mz.\\ (l + a(||z_||))(l-]l^.^(z)) 

< C'e-^^^'d^,{z,C+J ~ C'e-^^''dm{z,C+) 

< e-^'*d{{m,z),C+). 

Here in the second line we used the fact that ||2;|| < C||2;_|| if the distance is non zero (with 
C depending on 62). In the fourth line we used the continuity of the cone field, m 1— )> C^. 

This proves (4.6). The last claim (4.7) is immediate from the construction of and 
the fact that E+ n = {0}. □ 

We now regularize d{p, C^)^ uniformly with respect to a parameter e. It will eventually 
be taken to be h/h, where h is a small constant independent of h. Lemma 4.2 and the 
arguments of [31, §4] and [35, §7] immediately give 

Lemma 4.3. There exists > such that for any t > t^, there exists a neighbourhood Vt 
of K^^ and a constant Cq > such that the following holds. 

For any small e > there exist functions j± e C°°{Vt U ipt{Vt)) such that for p e 

l±{p) ^ d{p,C±f + e, 7±(p)>e, 
± (7±(P) - 7±(<^t(p)) + Coe ~ 7±(p) , 
^^■^^^ 9«7±(P)=C^(7±(P)^-I"I/^), 

7+(p)+7-(p)~^^(p,i^T + e- 

Following [31, §4] and [35, §7] again this gives us an escape function for a small neigh- 
bourhood of the trapped set. We record this in 

Proposition 4.4. Let 7± he the functions given in Lemma 4-3. For some M 3> 1 inde- 
pendent of e, define 

(4.13) Go = log(Me + 7_) - log(Me + 7+) 

on a neighbourhood V of the trapped set K^^ . 

For any to large enough we can find a neighbourhood of Ui <^V of K^^ and Ci > 0, such 
that 

Go = 0(log(l/6)), d^Go = 0(min(7+, 7.)-^) = 0(6"^) , |a| > 1 , 
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and such that for p & Ui r\p''^{[—5,5]), 

(4 14) d^iGoiVtM) - Goip)) = 0(min(7+,7-)-^) = (^(e"^) , |«| > , 

d{p, K'f > Cie ^ Go{,ipto{p)) - Go{p) > 1/Ci . 

In the intermediate region between Ui and {x : w{x) > 0} we need an escape function 
similar to the one constructed in [14, §4] and [24, Appendix]. We work here under the 
general assumptions of §1.3 and present a slightly modified argument. 

Lemma 4.5. Suppose that X is a compact smooth manifold, p G S™'(T*X;'R) , w e 
5'''(T*X; [0, oo)), k < m, and that (1.9) holds. For any open neighbourhood Vi of K^^ , 
there exists eo > and a function Gi e 25, 25))) such that 

G{p) = for p in some neighbourhood of K^^ , 

(4.15) HpG^{p) > for Pi w-\{e^, oo)), 
HpGiip) > forpe p-\[-5, 5]) \ U w-\{eo, oo))) . 

Proof. Call Uq = w'^^if), oo)) and suppose p G p-^([-25, 2S\) \ {Vi U Uq). We first claim 
that there exist r± = r±(p), r_ < < r+, such that 

(4.16) Vt+{p) e C/o or </7t_(p) e C^o , 

(4.17) (^T±(p) e'l^iUC/o. 

(Here and below we use the notation ipA{p) = Wtip) '■ t £ ^}-) 

To justify these claims we first note that since p ^ K'^^, '^b.{p) HUq which implies 
that 

(4.18) 3T,, iprM^Uo. 

Assuming that Ti < we want to show that (pT2 (p) G U C/q for some T2 > 0. Suppose 
that this is not true, that is 

(4.19) (^(o,oo)(p)n(yiuc/o) = 0. 

Then for any — >■ 00, 

Pi = ftM e P~'([-25, 26]) \ (\/i U Uo), moo){Pj) n {Vi U Uo) = 0. 

By (1.9) the set p~^{[—26, 26])\{ViUUq) is compact and hence, by passing to a subsequence, 
we can assume that p j — )■ p ^ Vi U Uq. We have (pt{pj) ^t{p), as j — )■ 00, uniformly for 
\t\ < T, and it follows that ^[o,oo){p) n (Vi U Uq) = 0. For t > -tj, 

iftip,) = vt+t,{p) C v'[o,oo)(p) C p'\[-25,25]) \ {Vi U t/o), 

which means that ipt{p) ^ViVJUq for t > —tj — )■ —00. We conclude that 

Mp) nViUUo = =^ (Pr{p) G K^^. 



DECAY OF CORRELATIONS 25 

This contradicts the property p ^ T4, and proves the existence of T2 > such that (fiTiip) ^ 
Vi U Uo. We call T_{p) = Ti, T+{p) = T2. 

In the case Ti in (4.18) is positive, a similar argument shows the existence of T2 < such 
that v^T2(p) ^ (^1 U Uo) 7^ 0. In this case we call T^{p) = T2, T+{p) = 7i- 

For each p G p^^{[—26, 26]) we can find an open hypersuface Tp, transversal to Hp at p, 
such that, if ipT±{p) G ^0, then for p' G Tp, 

(fT± {p') e C/o, (fT^ {p') e 14 u C/o . 

Notice that the closure of the tube Vlp =^ v^(T_,T+)(rp) does not intersect K^^ . Using this 
tube, we construct a local escape functions Qp G C^{Vtp), with the following properties: for 
some Ep > 0, and an slightly smaller tube fl'p C ftp containing </?(r_,T+)(p), 

(4.20) Hpgp{p')>0, p' ^w-\{ep,^)), H^Qpip') > 0, p' G \ (^-^((e,, 00)) U V^). 

Here €p is chosen so that if (pT±{p) ^ Uq then <PT±{Tp) C w~^{{2€p, 00)). 

To construct Qp we take (t, m) G {T_,T^) local coordinates: {t,m) ^ ^Pt{iT^) ^ ^p- 

Suppose that (^t-{p) ^ f^Oj and that V'(T_,T_+7)(rp) C w~-'^((ep, 00)) and 9?(r+-7,T+)(rp) C 
Vi U ^/q. Choose Xp G C^((T_,T+)) which is strictly increasing on (T_ + 7,T+ — 7) and 
non-decreasing on (T+ — 7,7+). Also, choose G C^{rp) with = 1- Then put 

gp{i~Pt{m)) Xp(t)'0p(m). Since iJ^fifp = Xp(t)'0p(m), (4.20) holds. A similar construction 
can be applied in the case where <Pt_(p) £ 14, <pr+(p) £ C^■ 
Prom the open cover 

p-\[-5, 5]) \ (VI U t/o) c U {^P : P e p-^([-5, 5]) \ (Vi U t/o)} , 

one may extract a finite subcover Uj=i ■ The closure of this cover does not intersect 
K^^ , so that the function Gi(p) Y^^=i 9pl{p) satisfies (4.15), for eo = min^ Cp.. □ 

We conclude this section with a global escape function which combines the ones in 
Proposition 4.4 and Lemma 4.5. The estimates will be needed to justify the quantization 
of the escape function in §6. The proof is an immediate adaptation of the proof of [31, 
Proposition 4.6] and is omitted. 

Proposition 4.6. Let V , Ui, Gq and to be as in Proposition 4-4) Wi be a neigh- 

bourhood ofK^^ such that Wi d Ui, Wi U ^PtoiWi) d V. 

Take X G C^(V) equal to 1 inWiUiptQ{Wi), and let Gi be the escape function constructed 
in Lemma 4.5 for Vi = Wi. Then for any P > 1, G G C^{T*X;R) defined by 



(4.21) 



G = xC3rGo + C4log(l/e) d 
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where C3 and C4 are sufficiently large, satisfies the following estimates 

|G(p)| <C6log(l/e), a"G' = 0(e-l«l/2), |«|>l, 
pew, ^ GivtM) - G{p) > -Cr , 
^' ^ peW,np-\[-S,S]), d{p,K'f>C,e=^G{vM)-G{p)>2T, 

p ep-\[-5,5])\{W,Uw-\{e,,cx^))) ^ G{ipto{p))-G{p)>Cslog{l/e) 
In addition we have 
(4.23) ^^^PG(P) 



for some constants Gg and N,. 

5. Analysis near the trapped set 

In this section we will analyse the cut-off propagator 

(5.1) X"exp(-ziP//i)x", 

where x"" = Op'^^ix), X ^ n Si and suppx C {p : d{p, K^) < R{h/h)^} for some R>1 

independent of h, h. We could take two different cut-offs on both sides, as long as they 
share the above properties. 

Our objective is to prove the following bound (announced in (2.4)): 

Proposition 5.1. For any e > and M > 0, there exist Co > 0, ho, and a function 
h !->■ /io(^) > 0; such that for < h < Hq and < h < ho{h), 

(5.2) llx'e-^*^/ VIU^^L^ < Co h-"*/-" exp (-|i(Ao - e)) , < i < Mlog 1/h , 
where Aq is given by (1.20). 

Since e"**^/'^ is unitary, the above bound is nontrivial only for < t < |^ log 1/h < t < 
Mlogl/h. 

5.1. Dcirboux coordinate charts. We start by setting up an adapted atlas of Darboux 
coordinate charts near K^, that is take a finite open cover 

and symplectomorphisms Kj : Uj Vj = neigh(0, M^'^). The standard symplectic coordi- 
nates on Vj then appear as a local symplectic coordinate frame on Uj. We may choose the 
coordinates such that they split into 

X={x,y), E={^,7j), y,r;eM*, x, ^ G M'^"'^* , 
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such that the symplectic submanifold fl Uj is identified with /C fl T^- C R^'^, where 

That is, (a;, ^) is a local coordinate frame on , while (y,?]) describes the transversal 
directions. 

We also assume that for each p & f] Uj, identified with some (x, ^,0, 0) e Vj, the 
subspace {(a;,y,^,0) , y e R'^*} is e-close to the transversal unstable space dK,j{E^), while 
the subspace {{xq, 0, ^o, rj) : V ^ M.^*} is e-close to the transversal stable space dK,j{E~). 

We want to describe the flow in the vicinity of K^, using these local coordinates. We 
choose a (large) time to > 0, and express the time-^o flow (ft^ : Uj^ — > C/jj in the local 
coordinate frames, through the maps 

(5.3) Kj^jQ — Kj^ o (f^^ o : Dj^jQ — >■ Aj^jg , 

where Dj^j^ C Vj^ is the departure set, while Aj^j^ C Vj^ is the arrival set. This is deflned 
when ^toiUjo) n Uj^ ^ and such a pair jijo for which this holds will be called physical 

Below we will also consider the maps k'j^j^ representing the time-n^o flow in the charts 

5.2. Splitting e~'^/'* into pieces. We want to use the fact that the propagator e^'*"^/'* 
is a Fourier integral operator on M associated with (fto- To make this remark precise, we 
will use a smooth partition of unity (tTj G C^{Uj, [0, 1])) such that each cut-off Tij is equal 

to unity near some Uj (s Uj, and the quantized cut-offs Hj =^ Op)^(7rj) satisfy the following 
quantum partition of unity: 
J 

(5.4) n = Uj U* = I microlocally in a neighbourhood of K . 

We will then split e~'*°^/'* into the local propagators 

(5.5) T^,,/= ^,e-''°''^%„ 

which can be represented by operators on L2(R'^) as follows. We define Fourier integral 
operators Uj : L'^{X) — )■ L^(]R'^) quantizing the coordinate changes Kj, and microlocally 
unitary in some subset of Vj x Uj containing ^^(suppTrj) x suppTTj, so that 

(5.6) Vi, n,n* = UjU*UjU* + o{K^) , 

The local propagators T^j^^^ are then represented by 

(5.7) 7;uo'=^%n*^e-^*°^/'^n,„z^;. 

Notice that for an unphysical pair jij^, Tj-^j^^ = 0{h'^)i.^^L^. For a physical pair jijo, Tj^j^ 
is an Fourier integral operator is associated with the local symplectomorphism Hj^j^. 
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Prom the unitarity of e"'*"^/'* we draw the following property of the operators Tj'j. 

Lemma 5.2. The operator-valued matrix T =^ (7ij)jj=i,...,j; acting on the space L^(M'^) 
with the Hubert norm \\u\\^ = 1 1 '"ill 1,2? satisfies 

\\T\\L2(jStdy^L^(j^dy — 1 + 0{h) . 



Proof. Prom (5.6), the action of Tj^j^ on L^(R'^) is (up to an error 0{h°°) 1,2^1,2) unitarily 
equivalent with the action of T^j^^^ on L'^{X). Hence, the action of T on L^(R'^)"^ is equivalent 
to the action of on L^(X)"^, where is the matrix of operators (5.5). 

To obtain the norm estimate we follow [2, Lemma 6.5], put % L^(X), U = e~**°^/'*, 
and define the row vector of cut-off operators C = (nj)j=i,...,j. The operator valued matrix 
can be written as T'" = C*{U Ij)C . Its operator norm on C{l-L'^) satisfies 

II^IIV) = ll(^')*^1lw = \\c*{u®ij)cc*{u*®ij)c\\ci:H^^ 
^\\C*{UU.U*®lj)C\\a:H^y 

Egorov's theorem (sec (3.6)) and [43, Theorem 13.13] imply that n"*^ == UTIU* is a positive 
semidefinite operator of norm 1 + 0{h), with symbol equal to 1 + 0{h) near K^, and its 
square root \/n^, as well as the product -\/n^n \/n^, have the same properties. Hence, 

rliV) = II ((Vm®^j)c)*(Vm®/j)ciU(«.) 

= II (Vm ij) c ((Tm ® ij) * Wan) 

= ||VmnVm|U(«) = i + c»(/i). □ 



5.3. Iterated propagator. In this section we explain how to use the Tj'j to study our 
cut-off propagator (5.1). 

First of all, Egorov's theorem (3.7) applied to T = HjU*, B2 = allows us to write 
(5.8) UjU*x''^xJl^j^* + 0{h'^U)L2^L-, i = l,...,J, 

where the symbol Xj = X° '^7^ ^ Si{T*MJ^). 

We start from a arbitrary normalized state u e L'^{X), and represent the part of u 
microlocalized near through the (column) vector of states 

u = (iij),=i,...,j, Uj = UjU* u, \\uf = = n^i) + ^ih"°)\Mh ■ 

j 
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The equations (5.7) and (5.8) show that 

j jl ,30 

Jl Jo 

= E K TnjoXlujo + 0{h^U)L.^L^ . 

Jl Jo 

Similarly, for n > 2 the propagator e~"''**°^/'* can be represented by iteratively applying the 
operator valued matrix T to the vector u. By inserting the identities (5.4), (5. 6) n times in 
the expression -^w^-tntoP/h^Wy^^ ^j^g following 

Lemma 5.3. For any n e N (independent of h), we have 

Jl^-intoP/h^v.^^ E %^I^jW.-l---rjiJoX7o%o+a(/^^/^^)L^-.L^ 

(5.9) ^ , 

= E ^^rMl [(r)1,„,o Xjo %0 + 0„(/i^^^)l^-.L^ , 
JniJO 

where the matrix of operators, T, was defined in Lemma 5.2. 

5.3.1. Inserting nested cut-offs. In this section we modify the Fourier integral operators 
Tjij, taking into account that in the above expression their products are multiplied by 
narrow cut-offs xj. 

By construction of Xji there exists i?o > (independent of h, h) such that for any index 
j the cut-off Xj £ Si is supported inside the microscopic cylinder 

(5.10) 5ro(V/^)V^ = {(^'1/'^'^) ■■ \y\^ 1^1 < Mh/W} C T*R^ 

Fix some Ri > 2Rq, and choose a function x° G C^(M^°'*, [0, 1]) equal to unity in the ball 
{|^|, 1^1 < Ri}, and supported in {\y\, \fi\ < 2i?i}. Normal hyperbolicity implies that there 
exists A > 2 such that the cylinders B, (see (5.10)) satisfy 

(5.11) S'i(fi2R) (E 5iJA, 

for all < i? < 1 and any physical pair j'j. 

We then define the familes of nested^ cut-offs {x"}neN) {x"}neN as follows: 

(5.12) Vn e N, x"(y,^) = X°(l/A-^ryA-"), 

(5.13) x"(x,y,e,r^) x"(y(V^)^/^ #A)^/') e ^i(r*R'^), 



^Below we use the notation X° ^ X for nested cut-offs, meaning tliat x" = 1 iiear supp(x). 
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uniformly in n: the smoothness of actually improves when n grows. Prom the assumption 
Ri > Rq we draw the nesting x° >- Xj for any j — 1, . . . , J. Furthermore, the property 
(5.11) implies that 

(5.14) for any physical pair j'j, x"'~^^ ° i^j'j ■ 
Prom these nesting properties and from (3.7) we easily obtain the following 
Lemma 5.4. For any j — 1, . . . , J we get 

(5.15) {xTx7 = xJ + o(hnL--.L- , xJ (xT = Xj + 0{}f^)i^._,L- . 

In addition, uniformly for all n & — 1, . . . , J, 

(5.16) Tj,,{xT = {x^'^TTMxT + OCh^L^^L^ ■ 

We will actually only use n smaller than Mlog l//i for some M > independent of h, h, 
so our cut-offs x" will all be localized in microscopic neighbourhoods of JC when h ^ 0. 
Furthermore, for such a logarithmic time the number of terms in the sum in the right hand 
side of (5.9) is bounded above by J" < h"^ for some > 0. As a result, taking into 
account the above cut-off insersions, this sum can be rewritten as 

(5.17) Ue-'''*°^/^x'"u 

= E ^^r.Ul T,„,„_, (x-^)"' ■ ■ ■ T,,,, ixYTn^o ixTxA + Oih^L^^L^ ■ 

In the next section we will carefully analyze the kernels of the operators Tjtj {x^)'^ ■ 

5.4. Structure of the local phase function. To analyze the Fourier integral operators 
we will examine the structure of the generating function for the symplectomorphism Kj^jg. 

We start by studying the transverse linearization d±K{p) of the for a point 

p e /C n Dj^jg. In our symplectic coordinate frames, this transverse map is represented by 
the symplectic matrix 5'jij(,(p) — S{p) e Sp(2dt,R) given by 

(5-18) '^"^^'W^)^"'- 

The linear symplectomorphism S{p) admits a quadratic generating function Qp{y^,y^,9'), 
where 9' G M'^* is an auxiliary variable: the graph of the map ^{y^-,rf) ^ ^{y^^V^) — 
S{p)^{y^ ,7]^) can be obtained by identifying the critical set 

Cq, = {{y\y',e') : de:Q{y\y',9') = 0} C R''\ 

This critical set is in bijcction with the graph of S{p) through the rules 

ri'^d,.Q,{y\y',e'), t]' = -dyoQ,{y\y',e'), {y\y\9') e Cq^ . 

More structure comes from taking the normal hyperbolicity into account. Recall that our 
coordinates are chosen so that that and E~ are e-close to {q = 0} and {y = 0}, 
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respectively. (Here we identified with their images under dnj - see §5.1.) This imphes 
the existence of a symplectic transformation R{p) e Sp(2dt,R) such that 

(5.19) R{p){{rj^O})^E;, R{p){{y ^ 0}) ^ E; , R{p)^I + 0{e). 
Since d±K{p) = S{p) maps E'^ to E'^^^y the matrix 

(5.20) S{p) ''^^ R{k{p))-' S{p) R{p), pelCn Dj,,, , 
is block-diagonal: 

(5.21) ^ (^(^) 

The normal hyperbolicity (1.17) implies that, provided to has been chosen large enough^, 
the matrix A(p) is expanding, uniformly with respect to p: 

(5.22) 3 i/ > 0, Vp e /C, ||A~'(p)|| <e-' <1. 

More precisely, for any small e > 0, if to is chosen large enough the coefficient v can be 
taken of the form z/ = to(Amin — £), where Amin > is the smallest positive transverse 
Lyapunov exponent of ipt near K^. 

Combining (5.19), (5.20) and (5.21) gives 

(5.23) ^{p)-[^ 0(e2A + '^A(p)-i)J ' 

This explicit form, more precisely the fact that the upper left block is invertible, allows 
to use a special type of quadratic generating function: 

Lemma 5.5. //to is chosen large enough, for each p the generating function Qp{y^,y^,0') 
can be chosen in the following form 

(5.24) Qp{y\y',9')^q,{y\9')-{y',9'). 

For any point {y^,y^,9') on the critical set Cq^, the auxiliary variable 6' is identified with 
rf of the corresponding phase space point. 

The function qp{y^, rf) can be written in terms of a symmetric matrix H{p): 

(5.25) qp{y\v') = \{{y\r,'l^H{p){y\rP)), H{p) = {^'^^ ^J^) , i/^^ invertible . 
The matrix S{p) is related to H{p) in the following way: 

(5.26) Sip) = ( -H2iH22^ 



^Recall that n represents ■ 
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Comparing with (5.23) we see that 

(5.27) iJ?2 = ^2i = A(p)-i + 0(eA(p)'i), H,, = 0{e), H^^ = 0{t), 

uniformly with respect to p. The quadratic phase function Qp will be relevant when we 
consider the metaplectic operator M{p) quantizing S{p) in §5.5.4 (see also (3.22)). 

From the study of the linearized flow in the transverse direction, we now consider the 
dynamics of 

(5.28) K = ^jiji) ' ^ jiji) ^ ^ ^ ^jijo ^ 

along the trapped set - see Fig. 2. When no confusion is likely to arise we use the notation 
D, and A, for the corresponding subsets of /C. There we have no assumptions on the flow, 
except for it being symplectic. 

Possibly after reflning the covers Uj, each map k can be generated by a nondegenerate 
phase function ip = iljj^jf^{x^,x^,6) defined in a neighbourhood of the origin in R'^"''* x 

^d-dt X Rfe^ where < A; < n — see §3.3. 

Since the Uj have been chosen small, this map — ?■ F;; can be assumed to be injective. 
Notice that the values of ip away from are irrelevant. 

We now want to extend ip into a generating function of the map k, at least in a small 
neighbourhood of JC. The intuitive idea is to "glue together" the generating function ip for 
K, with the quadratic generating functions Qp for the transverse dynamics d±K{p). 

Let us consider the following Ansatz for a generating function ^ for k: 

(5.29) *(x\ 9; y\ y°, d') = ^{x\x'', 9) + 5^{x\ 9; y\ 9'), 

with an extra auxihary variable 9' e R'^*. To simphfy the notations we spht the variables 
into longitudinal and transversal ones: 

(5.30) p\\^{x\x\9), p^^{y\y\9'). 

Lemma 5.6. The transversal correction, 5'^{p\\, pi_) , defined in (5.29) satisfies 

5^{p\\,p^)=Qp^^{p^) + 0{{y\e'f), 

where Qp^^{») is a quadratic form of the same type as (5.24,5.25), which depends smoothly 
on p\\. If p\\ e corresponds to the point (p^; p°) G F;^, then Qp^^ = QpO. 

In other words, the quadratic forms Qp^^ extend the forms Qp to a neighbourhood of C^, 

Proof. Since K, is preserved by k and carries the map k, we may assume that for any 

the function 5'${p\\,») has no linear part in the variables p±. At each point p\\ e 
(identified with some p^ & K.), the quadratic part Qp^^{p±) generates the linear transverse 
deviation from k near the point p°, namely d_j_K{p^). This means that = QpO, which 
has the form (5.24). This form could be attributed to the geometric fact that the graph 
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of d_i_K,{p^) admits {y^,rf) as coordinates (that is, this graph projects bijectively onto the 
(yS?7°)-plane). 

This projection property locally extends to the graph of k: in some neighbourhood of /C, 
the points of can be represented by the coordinates {p^ = (x^,,^") G )C]y^,rj^), where 
y^,r]^ G neigh(O). This property shows that 5^ can be written in the form 

(5.31) 5^if(p\i,p^)^S^p\^,y\e')-{y',e'). 

As explained above, the quadratic part ?p||(») of 5^(p||;«) must be equal, for G C^, to 

the corresponding QpO generating S^p'^). The equations for Cq, show that, if we fix small 
values {y^,9' — rf), then value py such that (p||, y^, |/°, € C"* is 0((|/^, 7?°)^)-close to 
C^. □ 



5.5. Structure of the propagators Tjij. Prom the above informations on the phase 
function ^ = ^j'ji we can the integral kernel of T = Tjij defined in (5.7) and quantizing 
the map Kjij, as an oscillatory integral. The general theory of Fourier integral operators 
(see §3.3) tells us that its kernel takes the form 

(5.32) r(x\y^x°,y°)= / X,,. «(P|h P±) e^^^^'i'^^^ + 0{K-),.^l. , 

where we use the notation (5.30). Let us group the variables {x, y) = X, (^, 77) = S, (6*, 9') = 
6. We may assume that the symbol a{X^, X^, 6) is supported in a small neighbourhood of 
the critical set C^. In particular, for small values of the transversal variables p±, a{»,p±) 
is supported near C^. From (5.7), this Fourier integral operator is microlocally subunitary 
in Vf X Vj. 

5.5.1. Using the cut-ojf near )C. We now take into account the cut-offs (x*^)^, and study 
the truncated propagator T {x^}"' appearing in (5.17). 

Lemma 5.7. For any k > we have 

(5.33) T {x'r = T^" + 0{h^)L^^L^ , 

k 

where the Schwartz kernel of the operator is given by 

(5.34) T^\x\y';x\y') f -^^^ a{p^^, p^) x'^'^'\y') xHv^v') e'^^^''^^"'-\ 

J {2'Kh) 2 

where x}^ =^ X^=oj given in (5.12). 

Proof. As in (5.16), the nesting property ;i^tt(fc+i) ^ y^k ^ implies that 

(5.35) T {x^Y = T {x^Y + 0(^~) , 
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uniformly for all A; > 0. The Fourier integral operator calculus in the class Si presented in 
Lemma 3.1 has the following consequence: 

which combined with (5.35) gives (5.33). □ 

5.5.2. Rescaling the transversal coordinates. Since we work at distances ~ (h/h)^ from the 
trapped set, it will be convenient to use the rescaled transversal variables 

(5.36) y^(h/hy^y, fj ^ (h/h)"^ri. 

Our cut-offs x*^; X*^ defined in (5.12,5.13) are then related by x'iVjV) = X'^ViV)- This 
change of variables induces the following unitary rescaling T : L'^{dxdy) LP'{dxdy): 

(5.37) Tu{x, y) = {h/h) ^ u{x, {h/h)^) = (h/h) ^ u{x, y) . 
We recall (see for intance [43, (4.7.16)]) that 

Ta'"{x,y,hD^,hDy)T* = a'"{x,y,hD^hDy), d{x,y,^,fj) = a{x,y,^,r)). 

k 

Through this rescahng, the operator is transformed into 

f x'= def ^j.x Y* : L^dxdy) ^ L^dxdy), 

with Schwartz kernel 

f-\x\y\x\f)= j j -^^j:f^a{p,,{h/h)h^) 



(5.38) W« (27r/i)^^ (27r/i)^* 

X X^^'^^^Kf) x''{f:V°) es'^^^ii^+^^^^ii 



5.5.3. Transversal linearization. The factor x'*('^+-'^)(y-'^)x''(y°, ?7°) appearing in the integrand 
(5.38) allows us to simplify the above kernel. Indeed, it imples that the variables p± = 
{y^^tp^rf ) are integrated over a set of diameter ~ RiA^. One can then Taylor expand the 
amplitude and phase function (5\E' in (5.38): 

Since we will restrict ourselves to values k < Mlogl/^|, uniformly bounded with respect 
to h, we may omit to indicate the /c-dependence in the remainder. As a result, up to a 
small error we may keep only the quadratic part of 5^, namely consider the operator with 

the Schwartz kernel 

de dfj^ 



(27r/.)'^ {27:hY^ 
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Combining the above pointwise estimates with the fact that a e S{T*'M?'^), and with (5.35) 
and Lemma 5.7, gives 

f x'= ^ Tix'^riy, hDy) + 0{K^)l^^l^ , 



uniformly for k < Mlog l/h\. 



5.5.4. Factoring out the transversal contribution. For each G suppa(«,0), the quadratic 
phase (5p||(*) generates a symplcctic transformation S'(p||) (which, in the case py G 
corresponds coincides with the transformation S'(p°) of (5.18)). As already shown in (3.22), 
this phase allows to represent the metaplectic operator M(p||) : L'^{dy) — )■ L'^{dy) which 
/i-quantizes this symplectomorphism: 

(5.40) M{p\\){y\y') "^'{27:11)-'^ [ det{H,^{p\\)Yl^ ei'^n^'^Uf,\ 

where Hi2{p\\) is the block matrix appearing in Qp^, similarly as in (5.24,5.25). 

Remcirk 5.8. In the expression (5.40) we implicitly chose a sign for the square root of 
det(i?i2(p||)). Indeed, the metaplectic representation of the symplectic group is l-to-2, a 
given symplectic matrix S being quantized into two possible operators ±M. The relations 
(5.27) and the uniform expansion property (5.22) show that det(ifi2(p||)) does not vanish 
on the support of the amplitude a(«, 0) (which is a small neighbourhood of x = 
y^ — ff — 0}), so we may fix the sign in each connected component of this support. This 
remark will be relevant in §5.6. 

Defining the symbol 

def a(P||,0) 



|2 ' 



det(i/i2(p||))^ 

we interpret the operator T in (5.39) as a Fourier integral operator with an operator valued 
symbol, M(p||), where M is given by (5.40). That fits exactly in the framework presented 
in Proposition 3.5: 

f{u ® ^;)(x^ y^) = (27r/i)-^^ f [ a{p\\) [M{p\\)v]{y^) e^'l'^P\\\{x^)dx^de. 
We now apply Lemma 3.4 to see that 

(5.41) f = Op)f (M)TII + Of^{h)T,mr.+e_^V^ , 

where the Schwartz kernel of is given by 

(5.42) rll(a;°,x^) = (27r/i)-*^ /" a{p\\) e^'^^'^iO dO , 



36 S. NONNENMACHER AND M. ZWORSKI 

while the operator valued symbol M{p^) is the metaplectic operator /i-quantizing S{p^), 
where — and S{p^) is given in (5.18). We summarize these findings in the following 

Proposition 5.9. Suppose that the Schwartz kernel of T is given by (5.32), x'') x'^ 
given in (5.12), and T is the unitary rescaling defined in (5.37). 

Then for k < K{h), where K{h) may depend on h hut not on h, 

(5.43) TiTix^'T) r* = Ovl{M)T\x''r + Oih^)^.^^. + Oj^ih)^^^^^ , 

where T" is given by (5.42) and M{x^,^^) given by (5.40) with py G determined by 
{x^,^^) = {x^,dxiilj{p\\)). Here and below we use the abbreviation (x^)"' == {x^)'^{y,hDy). 

Proof. Lemma 5.7, (5.39), and (5.41) give (5.43) with the remainder 

0(11°°) L2(dxdy)^L2(dxdy) + O -^{k) l2 (^dx)®V^^ ^ L'^ {dxdy){X^)^ ■ 

The definition of x^ in (5.12) and (3.24) show that 

and that gives the remainder in (5.43). □ 

5.6. Back to the iterated propagator. We can now come back to (5.9) and (5.17), 
re-estabhshing the subscripts jk+ijk on the releavant objects. We rescale all the operators 
by conjugating them through T. Fixing the limit indices jo,jn, we want to study the sum 
of operators obtained by conjugation of terms in (5.17) by T: 







(5.44) 



r[T-h.joixTr* = T[J2 n TM{xT]r* + oChnL2^L2 



j k=n—l 



where the sum runs over all possible sequences j — ■ ■ -ji- A sequence (which could 
be thought of geometrically as a path) jnjjo will be relevant only if it is physical, meaning 
that there exists points p G such that (fktoip) ^ ^ife ^oi all times k = 0, . . . ,n (we say 
that the path jnjjo contains the trajectory of p). Any unphysical sequence leads to a term 
of order 0{h°°). On the other hand, for a given point p G there are usually many 
sequences j containing its trajectory, since the neighbourhoods (C/j)'s overlap, and so do 
the cut-offs (tTj). 

We now simplify the expression (5.44), in the following way. 
Lemma 5.10. In the notation of (5.9) an(i(5.44), and for n < Mlogl/h, 
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Here T^jj^ is a Fourier integral operator on L?'{dx) quantizing the map : Vj^ — >■ Vj^, 
defined on the departure set Dj'^j^. For each p e ^j„jQ ~ '^"(-^jnjo) (^^^ arrival set) the 
operator valued symbol -^j^jo(p) is a metaplectic operator quantizing the symplectic map 

Sl,,{{nn-\p)) = d^K-{{R-)-\p)). 

Proof. If wc insert the approximate factorizations (5.43) in a term j of the sum in the left 
hand side of (5.44), this term becomes 

(5.46) OpnM,„,„_j7;ll ,_^(X"-^)'^ • • • OvmnjoKiM'T + 0(^~)l^^l- 

We now observe that just as we inserted the cut-offs to obtain (5.17) from (5.9) we can 
remove them so that each term becomes 

(5.47) OpnM,„,„_ J7;|l • • • 0^l{M,,,,)Tl^{x'r + (^(^~)l^-.l- 
We can now apply Lemmas 3.3,3.4 and Proposition 3.5 to see that 

where we use the shorthands 



(5.48) 







def 


rP . 


T"!! . . . T"!! 


JnJn—l' 


■■jo 




Jnjn—l 


jn-ljn-2 jlj'o' 


Mi 

Jnjn—l 


■■jo 


def 


^'^^^jnjn- 


-l)(-^?n-lJn-2 ° ^jn-ljn) 






def 








■■jo 




'^jkjk-1 


° ^jk-ijk-2 • • • ° ^jija- 



These expressions only make sense for physical sequences jnjjo- The map Hj^jjo is defined 
on the departure set Dj^jj^. 

The metaplectic operator Mj^jj^{p) quantizes the symplectomorphism Sj^jj^{p'^), with 
p — 'k^{p^) e ^j„jjo- This symplectomorphism represents, in the charts Vj^ — )■ the 
transverse linearization of the flow ipnto at the point K^^^{p^). As a consequence, the sym- 
plectic matrix Sj^jj^lp^) is identical for all sequences jnjjo containing the trajectory of p", 
and we call this matrix S^^j^{p'^). Hence, two metaplectic operators Mj^jj^^p), Mj^j'j^{p) 
corresponding to two different allowed sequences can at most differ by a global sign. 

For all p in the arrival set 

^ho = U ^™ = '^:>n({p e n K\ cp-ntoip) e U^o}) , 
j 

we choose the sign of the metaplectic operator Mjl^j^{p) quantizing S^^j^{p^), such that 
Mj^^-^(p) depends smoothly on p on each connected component of A^^^^ (there is no ob- 
struction to this fact, due to the property mentioned in the Remark 5.8: the symplecto- 
morphisms S^^j^{p) also have the form (5.23)). Hence, for each physical sequence jnjjo we 
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'3231 jo 




''323'-l3o 



Figure 2. Schematic representation of the departure and arrival sets for 
j of length 1. We show two physical sequences j2jijo and j2jijo and the 
corresponding maps (5.28). As remarked there we use the same notation for 
the departure and arrival sets on /C. 



have 

(5.49) ^™(p)=£™(p)^Lo(P)' P^Dj^iJo, 

for some sign Sj^jjo{p) £ {±} constant on each connected component of Aj^jj^. As before, 

the functions p i— )■ £j„jjo{p), p ^ ^jnjjoip) "^an be smoothly extended outside Aj^jj^^, into 
compactly supported symbols. Lemma 3.3 and the identity (5.49) give 

(5.50) OpnM,„,„_,...,jr,ll ..^ = Op-{Ml.J (£,„,-,J- tI-.^ + 0~,{h),.^,^^^^r.,^_,,. . 

When (x°)^ is inserted in (5.48) and (5.50) we apply (3.24) to see that 

0(/i-2"--<^*/i)^,(^^^^^.„.,_,^^^,(x°)* = 0(/i-^~<'*/i)^.^^. = 0~,ih)r..^^., 

and hence that error term can be absorbed into 0{h°°). 
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Returning to (5.46) we see that the sum in the right hand side of (5.44) can be factorized 
in the following way: 



(5.51) 



'0\w 



ovmi,:) ( (.,„,,,)-) {^r + o{h-),.^,. , 



with a uniform remainder for n < MIoet 1/h. Let us put T"",- =^ • Tf ■■ (e,- oo,,)™, so that 
the above identity reads exactly like in the statement of the Lemma. The operator T'^^-^ 
is sum of Fourier integral operators T^^^-^ defined with different phase functions ipj^jjo , yet 
these phases generate (on different parts of phase space) the same map : D^j^j^ — > ^l«jo ■ 
Hence, T^j^^ is an Fourier integral operator quantizing k". This completes the proof of 
(5.45). " □ 

The next lemma shows that the Fourier integral operator T^j^j^ is essentially subunitary. 

Lemma 5.11. Let M > 0. For any small h > 0, there exists ho — ho{h) such that, for 
any sequence j of length n < Mlog 1/h and any h < ho(h), the operator T^j^^ satisfies the 
following norm estimate: 

(5.52) \\T^!jLHd.)^L^d.)<l + 0{h). 

Proof. We first note that we can bound the left hand side of (5.52) by h~'"^ , for some C 
— that follows from a trivial estimate of the terms t]' ■ • in (5.51). 

To prove (5.52) it is clearly enough to prove the bound ||Tj^'io(x°)'"||L2(dxdif)^L2(da;djj) < 
\^0{}f^). From Lemma 5.2 we know that ||T"||(^2)j^(£,2)j < l + C)(/i), which implies that 
ll[T"]jojnlU2^i2 < 1 + 0{h). Lemma 5.10 then shows that 

(5.53) l|OpnMjl,o)7;l(x°)1U^-.^. < 1 + 0(h-^). 

The family of unitary metaplectic operators p h- )■ Mj^jy(p)~^ is well defined for p in the 

neighbourhood of the arrival set , and T"",- is microlocalized in any small neighbour- 
hood of A" . X , C Vi„ X Vi„. Lemma 3.3 and (3.24) then show that 

Jnj{] J 71 JO ^ / 

T;!,,{^r - opn(Mj:^.j-)opnMj:^.j7;i(x°)^ + omT;^j)LHa.)^v^-^-^.-.LA^r 

where we used the above a priori bound on ilTf'Lll- 
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Just as before we can insert the cut-off x" (see (5.12)) witfi a 0(h°°) loss. We also 
introduce a cut-off ip — ip{x,^) to a small neighourhood of Aj^j^. (It was not necessary 
before as Tjj^j^ provided the needed localization.) This and and (5.53) give the bound 

Since by Lemma 3.3 and (3.24) 

[opn(^i:.o)"V)(x")i*opn(Mj:,,j-V)(x")^ = [rrrnxTTixT + o-,{hu^-,L^, 

we have 

l|Opn(^Lo)"V)(xTII < \\r\\\\{xT\\ + o-,{h) <i + o{h), 

and the bound (5.52) follows. □ 

5.7. Inserting the final cut-off. We now return to the operator -i^"'e~**"°^/'^x"'. From 
Lemma 5.3 we easily obtain 



(5.54) 



X~ e T^-^u = > II,///: vr III r"L- vr + tvi A),2/7,2 

jn ,30 



where in the first line we used (5.8), while in the second hne we used (5.15). Hence our last 
step will consist in estimating the norm of the operator (x*^)^ [T""] j„jo (x°)"' (or its conjugate 
through T). To this aim we will use Lemma 3.4, Proposition 3.5 and the factorization (5.45) 
to obtain 

Here the operator valued symbol -/V"j^jo(p) = -^j^jo(('^")~^('^))' P ^ ^7njo^ ^ metaplectic 
operator quantizing the symplectic map S'j^j^{p) — d±K'"{p). (Having it on the right now 
makes the notation slightly less cumbersome.) 

In Lemma 5.11 we control the norm of T"''- . There remains to control the norm of the 

JnJQ 

factor (x^)"" ^Phi^j" jo) (x°)'"- For that it is enough to control the operator-valued symbol 

5.7.1. Controlling the symbol. In (5.19) we defined, for each point p & K. D -Djij,,, a sym- 
plectic transformation R{p) G Sp{2dt, M) which maps the y-spacc to E'^ and the ?7-space to 
. This transformation is e-close to the identity and in particular it is uniformly bounded 
with respect to p. 
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By iteration of this property, for any po £ -D? ' 



is block-diagonal in the basis {y,r)): 

(5.56) = ('^"f"") «^,_.) , 

where A"(po) is expanding. We may quantize R{p) into metaplectic operators A{p), and 
define _ 

NLoiPo)-MPnr'Nl,^{po)A{po) 

which quantizes S'j^j^{po). 
We can then rewrite 

We are interested in the — > norm of this operator. Since metaplectic operators are 
unitary, and using the covariance of the Weyl quantization with respect to metaplectic 
operators, this norm is equal to that of 

i^pT iZo ° si,,{po)-'r , = ^ ° Ripn), = x° ° Ripo)- 

The block diagonal form of S'j^joiPo) shows that 

(5-57) o {SlJpoT'Uv) = X°o(A"(po)-^^/A"(po)r/) . 

We may now invoke the following simple 

Lemma 5.12. Suppose that A is amxm real invertible matrix and that xi, X2 £ ^(M^"*). 
Then 

(5.58) \\xr{x,hDx)X2{AxjA-'hD,)\\L2^j^rn)^L^^rn) < C|detA|U-^, 

where C depends on certain seminorms of xi o-nd X2, but not on A. 

We remark that the upper bound becomes nontrivial only if |detA| <^ /i™/^. In that 
case, one cannot apply the /i-symbol calculus any longer (which would predict that the 
norm is equal to max^. ,? |xi(a;,^) X2{Ax A~^ ^) \ + 0{h)) because the second factor is not 
a symbol in the class 5'(M^"*). 

Proof. If we put Xj{x,Z) ^ Jj^^ Xj(x, ^)e*^'^'^^d^, then the kernel of the operator in the 
lemma is given by 

K{x, y) = —\— [ xi e) X2 [^M-^v] e^(--'«)A+^(-^'''>/^de drj dz 
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We will estimate the norm using Schur's Lemma and hence we need to show that 
(5.59) (^max J \K{x,y)\dy^ ^max J \K{x,y)\dx^ < C^l det^l 

Making a change of variables Z = {x — z)/h and X = {x + z)/h we obtain 

\K{x,y)\dx <Ci{meix\x2\)\detA\h~'^ II \xi{X, Z)\dZdX < C\det A\h-'^ . 



I 



To estimate the integral in y let 

F{Z) = max |xi(., Z)l G{Y) = max |x2(«, Y)l 

noting that our assumptions give F{Z) = G{Y) = 0{{Y)-°^). Changing 

variables to Z — {x — z)/h and Y — {Az — Ay)/h we obtain, 

j \K{x,y)\dy ^CsJJ F{Z)G{Y)dZdY < C. 

This proves the upper bound (5.58). □ 

Applying Lemma 5.12 to the product (5.57) we get the bound 

Since the transformations R{p) are uniformly bounded, the prefactor C(XpQ,Xp„) is uni- 
formly bounded with respect to po- On the other hand, the determinant of A"(po)~^ can 
be bounded as follows. 

Lemma 5.13. Take £ > arbitrary small. Then there exists Cg > such that, 

Vn > 1, Vpo e Dl^^, I det A"(po)-'| < C,e-(^°-^)"*» , 
where Aq was defined by (1.19). 

Proof. This follows from writing the definition of Aq using the local coordinate frames. □ 

We have thus obtained the following upper bound: 
(5.60) \\{)^f Nl^M) imLHrn-^LHay) < a ^"'^^/^ e-(^°-)"*° , 

valid for any n > 1 and any po £ -^j„jo' particular, the time n may arbitrarily depend 
on h. 

When n < Mlog 1/h, for M > arbitrary large but independent of h or h, we combine 
this bound with (5.17), Lemma (5.11) and Lemma 5.10 to obtain the estimate (5.2), which 
was the goal of this section. 
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6. MiCROLOCAL WEIGHTS AND ESTIMATES AWAY FROM THE TRAPPED SET 

In this section we will justify the estimates described as Step 2 of the proof in §2. 
That will involve a quantization of the escape function G given in Proposition 4.6 with 
e^{h/h)2. That means that we will use the calculus described in §3.2. 

6.1. Exponential weights. Suppose that g e C^(T*X;M) satisfies the following esti- 
mates: 

(6.1) ^^^^^<C{l + Ch/hy^d{p,p')f, d;g = 0{{h/h)-\-\l'), H>0, 

9\P ) 

for some and C , and for some distance function d{p, p') on T*X x T*X (since g is 
compactly supported, the estimate is independent of the choice of d). We note that G 
defined in Proposition 4.6 with e = {h/h)^ satisfies these assumptions. 

We first recall a variant of the Bony-Chemin theorem [5, Theoreme 6. 4], [43, Theorem 8.6] 
in the form presented in [31, Proposition 3.5, (3.21), (3.22)] (as usual g'^ — Op^{g)): 

Proposition 6.1. Suppose that g G C^{T*X) satisfies (6.1). Then 

(6.2) exp(^7-) = 6-, 
where the symbol b{x, ^) satisfies the bounds 

(6.3) \d%{p) I < Ca e^^P^ (h/h) , 
in any local coordinates near the support of g. 

Ifsuppg U, for an open U (s T*X, then 

(6.4) d^d^^{b{x,0-i)^O{h^{0-n, 

Also, ifAe *^°'"P(X), B e ^^"''(X) and C e *?''"P(X) then 



2 

r , 1 



e'' ^ A + i{hh)-^Ai, Aie^i{X), WFh{Ai) cWFh{A), 

(6.5) e^'"Se-^"' = S + Mi, B.e^'r'^iX), WFh{B^) cWF^iB), 
gP^C'e-^"' = C + iUCu Ci e ^r"^P(X), WF,,(Ci) c WF^C). 

2 

The assumptions in (6.1) show that expg^ is an order function for the 5*1 calculus - see 
[31, §3.3, (3. 17), (3. 18)]. Hence we can apply composition formulae. In particular if gj, 
j = 1,2 satisfy (6.1) then 

(6.6) exp(^r) eM92) = |a"c(p)| < G^expig, + g^){h/h)-^''^'\ 

Because of the compact supports of g^^'s and because of (6.3) derivatives can be taken in 
any local coordinates. 

The consequence of (6.6) useful to us here is given in the following Lemma. 
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Lemma 6.2. Suppose that A e ^?'"^(^) and that 

2 

a{A)^a + 0{{hhy^)^^, aeC^{T*X)r\Si{T*X). 
IfUf^-^ = {p e T*X : (i(p,suppa) < {h/h)^}, then 

(6.7) \\Ae^^ e^^ = sup{\a\e^'+^') + OCh sup e^'+^^) + O{hhog{l/h)). 

Proof. We first consider this statement in M"^. We apply the standard rescahng (3.4) noting 
that (6.1) imply that fhj = Qj are order functions. If d is the Euclidean distance and if we 
put 

UNiP) = (1 + dip, t/))-^, U (h/hy^U,;,, 

then un is an order function for any A^, and d G S^un) for all A^. We have 

A — Op'^{a+ {h/h)^ai), for some oi e 
and hence, after rescaling, 

ieOPl'(^i)e°Ps'(52) ^ Op^(6) + {h/hy^Opl{h), 
b e S{nN'rhirh2), b — ae^^"*"^^ e hS{nNThim2), bi e 5'(?fiim2). 

Put 

M = M{h, h) = supnArmirfiz < sup ((l + d{p, U))~^ e^^^^^^^'^^^A 

< (^supe^^+^2) (^l + sup(l + CiC2rf(p,^))-^+^i+^2^ 

< Csupe^^+^^ 

where we took N > Ni + N2, with Nj, Cj appearing in (6.1) for gj. 

We now apply [43, Theorem 13.13] (with h replaced by ^) to 6/M e 5". That gives 

||0p;^(6)|| = sup |a|e»i+^2 + 0(/i) sup e^i+»^ 

Since fhifh^ — O(log(l//i)), applying the same argument to 6i/log(l//i) gives (6.7). 

The calculus is invariant modulo 0{{hh)^) terms (sec (3.5) and [11, §5.1], [42, §3.2]), so 
these local estimates on imply similar estimates on manifolds. □ 

The next result is a version of (3.6) for exponentiated weights g. It is a special case of 
[31, Proposition 3.14] which follows from globalization of the local result [31, Proposition 
3.11]. We state it using concepts recalled in §3.3. 
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Proposition 6.3. Suppose that T e I^'^p{X x X, T'^) where k : C/i ->■ U2, Uj C T*X, is a 

symplectomorphism, that g e C^{T*X) satisfies (6.1), and that A e ^^°™p. Then 

2 

^^■^^ neir'^iXxXX.), B,Ce^.{X), a{B) = K*a{A). 

6.2. Estimates away from the trapped set. We now provide precise versions of the 
estimates (2.5) and (2.6) described in the Step 2 of the proof in §2. 

For the escape function G constructed in Proposition 4.6 we define the operator 

(6.9) Op^G) e \og{h/h)^r'^^{X), a{G) = G + 0((M)^-). ■ 

2 

Since G satisfies (6.1), Proposition 6.1 describes the exponentiated operator e*-^™ = e'^'^^i^). 
The first lemma shows that the weights are bounded near the trapped set: 

Lemma 6.4. Suppose that x e C^{T*X) D Si{T*X) has the propeHy 

(6.10) supp X C {p e T*X : d{p, K^^) < Go{h/h) ^ }, 
for some constant Go > 0. 

Then for some constants ho,ho,Gi > we have for < h < Hq, < h < Kq, 

(6.11) llx^'e^'^'ll < Ci, ||e^"x"|| <C'i. 

Proof Since ?(x'^) = x + 0(^^^^)5i , and \G{p)\ < G2 for d{p,K^^) < (Go + l)(VM^ the 
estimates in (6.11) follow directly from Lemma 6.2. □ 

The main result of this section provides bounds for the conjugated propagator. It relies 
heavily on the material about the propagator for the complex absorbing potential (CAP) 
modified Hamiltonian, exp{—it{P — iW)/h), presented in the Appendix. 

Proposition 6.5. Suppose that G"' is given by (6.9) and that A e ^^°°^p(X) satisfy 

(6.12) WF,(A) c p-\{-S, S)) n w-\[0, eo)), 
for some eo > 0. 

Then for some constants hQ,ho,Gi > we have for < h < Hq, < h < ho, 

(6.13) ||e-«"'e-^*°(^-^^)/^«"'A|| < e^^^ 
If X satisfies (6.10) and in addition 

(6.14) x(p) = 1 ford{p,K^') < \Co{h/~h)h, \p(p)\ < 5, 
where Cq is a large constant dependending on to, then, if \\A\\ < 1, 

(6.15) 11(1 - x")e-^"e-^*«(^-'^)/'^e«'"^|| < e'^ 
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where F is the constant appearing in the definition (4.21) of G 



Proof. Let A^g e'~^'^ Ae Then (6.5) in Proposition 6.1 shows that 

(6.16) A_G^A + OL^_,L.{h^) = ^(l)^^^^^ and A_g = AA_g + 0{h'=°), 

where A satisfies (A. 10). To prove (6.13) we use the notation of Proposition A. 3, and 
rewrite the operator on the right hand side as 

^-G-^-ito{P-iW)/h^G-j^ = e-«"'e-^*°^/'^e«""e-«>;j(to)A_Ge^"' + 0{h^L.^L^ 
^^■^^^ = e-«"e-*«^/^^"'C(to) + 0(/i^)L^^L^ , 

where using (6.5) and Proposition A. 3, 

C{t^) e ^T^iX), W¥h{C{t^)) c WF^(^) n w-\Q). 

2 

Since 

e±«'" = Be^^" + {I-B) + 0{h'^)L2^L^, for some B G ^'=°'"p(X), 
Proposition 6.3 (apphed with A = I) and (3.8) show that for some Bq e ^i(X), 

From this and (6.17) we see that to prove (6.13) it is enough to show that 

e-KG)"" e^^Bi = 0{1)l.^j^., B, e vl/--P(X), 

WF,(Si) C p-i((-5, (5)) n ^-^([0, eo)). 

Lemma 6.2 apphed with gi = —fl^G and g2 = G, and the property G — <Ptf^G < Cr in 
(4.22) which holds in a neighbourhood of WFh{Bi), give (6.18) and hence (6.13). 

To obtain (6.15) we proceed similarly but applying the property ip^^G — G > 2T which 
is valid outside a (h/h)^ neighbourhood of - see (4.22). In more detail, Proposition 6.3 
applied with A = 1 — gives 

(1 - x'")e-^"'e-^*o^/'^e«'"e-^>^(io)A_Ge^™ + 0{h'^U._,,^ = 
^-UoP/H,-i.;,Gr _ (^;^^)-)e-^-Vj(to)A_Ge^- + 0{h^,.^,. , 

where we used the boundedness established in (6.13) to control the lower order terms. 
Defining xi ftoX^ we have, by the invariance of under the fiow, 

XI = 1 for dip,K')<C^ih/h)K b(p)| < 5. 
Let ip e C^{T*X) be equal to 1 in the set Wi of Proposition 4.6, and suppV' C (w"^(0))°. 
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Since (6.5) and Proposition A. 3 give 

||e-«>^(to)A_Ge«ni < l|e-^>l(^o)e^niP|| < \\A\\ {\\A\\ + 0,.^,.iV')) 

it is enough to show that 

(6.19) ||e-KG)'"e^'"(l - xr)r\\ < e-'^/^ 

(6.20) ||e-K/')'"f/^"'(l - xi )i^ - i'")B\\ < Ch, 

for B e ^™°ip(X) with WF^iB) C W2 where W2 C w'\[0, eo/2]) f] p-\[-5, S]), is as in 
Proposition 4.6. Both inequahties foUow from Lemma 6.2 and properties of G in (4.22). 
For (6.19) we apply (6.7). For (6.20) we note that on 

vlG-G>CslogCh/h), on supp(l - ^) n 1^2, 
and (6.7) gives the estimate. □ 



7. Proof of Theorem 2 
We first prove (2.2) which we rewrite as follows 

(7.1) ||C/S^||l2(x)^l2(x) < Ce-"*°(^°-^)/^ M.log i < n < Mlog i 

h h 

where 

with to chosen in previous sections, and 

(7.2) ^e*^°'"P(X), WFh{A) Cp-\{-6,6)). 

To apply the estimates of the last two sections we first observe that Proposition A. 2 
implies that for any r there exist Bj e ^"^"^p, j = 1, • • • r, each satisfying (7.2), such that 

r 

(7.3) C/£A = J]C/gS, + 0(Ol2^l^ Br^A, 

where the constants in the norm estimate 0{hP°) depend on r. This means that, for r 
independent of h but depending on /i, UqA can be replaced by the product of operators 
UcBj, to which estimates of the previous section are applicable. 

We now want to decompose Uq in such a way that the estimates obtained in §§5,6 can 
be used. For that we define 

(7.4) Ug^Ug,+ + Ug,-, Ug,+ ^^Ugx'", Ug,- = Ug{1 - xT ■ 
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We note that Proposition 6.1 shows that 



(7.5) 



where xf == X^^ ''^^^ w)/h^itp/h ^ Proposition A. 3 apphed with with P replaced 

bv -P. A e ^^°'"P satisfying WF;,(/ - A) n WFhix"") = 0. In the notation of (A.12), 
Xf = X"'VA{ty, VA{t)* e ^pP(X). From (A.12) 



a{VA{t)) = exp {-^-j\*_^W^ a{A), 



with a full expansion of the symbol in any coordinate chart given in Lemma A. 4. For p e 

suppx, d{p,K^) — 0{h^), and as is invariant under the flow d{(p_sip), K^) = Os{h^). 
But that means that on the support X) V^-s^ = for s < t, where t is independent of h, 
as long as h is small enough. This means that WF/i(/ — V^(t)*) fl WF/j(x) = and and 

hence, Xt = X + C(/i^)5i- 

Returning to (7.5) this means that for t < Clog(l//i) (in fact for any time bounded 
independently of h) , we have 

Using the notation (7.4) 

U^=J2^0,er.---UG,e,UG,e, 
eGS(n) 

where we used the symbolic words e = ei ■ ■ ■ G S(n) = (±)". Now, for each word 

e ^ • • • , call niie) ( nji{e), respectively), the number of consecutive (— ) starting 

from the left (the right, respectively): 
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Given integers nL,nR, call E(n, 71^,71^?) the set of words e e such that ^^(e) = ul 

and nij(e) = rii,. The decomposition (7.7) can be split into 

where the sum runs over ul, nR> such that 1 < til + riR < n — 1. 
We make the following observations: 

riL, nn) = {(-)"^ + (-)"«} , if + = n - 1, 

riL, ur) = {(-)"^ +e'+ (-)"« : e' e E(n - - riij - 2)} , iinL + nR<n-l. 
Hence, the above sum can be recast into 

n-l 

US = us,- + E (uc-r ug,-, {UG,-r''''-' 

(7.8) "i=o 

+ (Uc-r ug,+ ([/g)"-"«-"--^ ug,+ {UG,-r'' 

where the last sum runs over ni.nR > such that 1 < Ul + ur < n — 2. 

The following lemma provides the estimate for terms in the first sum on the right hand 
side of (7.8): 

Lemma 7.1. For h > h > small enough, the following bound holds for ro < r < 
Colog(l//i), r eN, 

(7.9) \\Ug,+ USUg,+ \\l,^l, < Ch-''^/' exp (-|tor(Ao - e)) , 
where the constant C is uniform with respect to h, h and r. 

Proof. Lemma 6.4 shows that e'^'^x^i X^^^^"" — 0{1)l2^l2. Also, Lemma 6.2 shows that 
for xi £ Si with the same properties as x but equal to 1 on the support of x, we have 

^.,^-G^ = X-e-^'-Xi + 0(h^)L2^L^, e^^x" = X^e'^'^x'' + (^(/^°^)l^^l- 
Using (7.6) the operator on the left hand side of (7.9) can be rewritten as 
Ug x"" {UgY^' x"" = Ug X- e-«™ ^-^i^r^^t^^l^ + 0{h^^.^^. 

= UG{x'"e-^'"Xi + e'(/i°°))e-'('"+')*°^/'*(xre^'"x'" + 0{}f^)) 

Hence, 

\\Ug,+USUg,+ \\l^^l^ < \\Ug\\ llx^e^^'ll \\Xie"^'^'''^'°''/'xr\\ He^'^x'H + 0(/^°°) 
< C \\xie~"^"^^^'''^^'"Xi\\ , 
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where we used the fact that the operators Uq, e*^"'x"' and x"'e~*^"' are uniformly bounded 
on L^. We can now apply Proposition 5.1, replacing i by (r + l)tQ and x by Xi- ^ 

Let us now take n = Colog 1/h, with Cq ^ 1. We recall that T in (2.5) was assumed to 
satisfy F > toXo/2. We will use the bounds (7.9), and Proposition 6.5: < 

Returning to the estimate for UqA we first observe that (7.3) and the estimates (6.15) 
in Proposition 6.5 give 

(7.10) \\Ug,-M < e-'"^ + a(/i^). 

In (7.8), for each £ = 1, . . . , n — 2, we group together terms with nL + riR — £, and apply 
Lemma 7.1 and (7.10): 

n-2 

> An 



n-2 _ 



i=i 

An-e 



By taking Cq — 1/e we may absorb the prefactor /i"*^*/^ and obtain, for /i > small 

enough, 

(7.11) \\U^A\\<C e^p{-nto^^^^), n^M,\ogl/h. 

We can now complete the proof of (1.18) following the outline in §2. We first note that 

(7.11) gives (2.2), so that (see (2.3)) for 

(7.12) z e [-5/2, 5/2] - ih[Q, (Ao - 3e)/2] 
and A e *^°™p(X) satisfying (7.2), 

(Pg - z)Qa{z) =A- R{z), R{z) = OCh)L^^L^, 

^^^^^ Ji e-'^^^^-^y^Adt = O , T{h) = M, log 1/h. 

We now apply this estimate with A e *^°'^p(X) such that a{A) = 1 in p-^ {-35 /A, 35/4:)p 
w-^([0,eo)). Then Pq - z e *T(X) is elliptic outside of WFh(^). Hence, using the *i 

2 ^ 2 

calculus of §3.2, there exists Qa{z) e 4fJ"^{X) such that 

2 

(Pg - z)Qa{z) = I-A + R{z), R{z) = 0{h)L2^L2. 
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The Predholm operator Pq — z has index since Pg + i is invertible for small h. It follows 
that for h small enough and z satisfying (7.12) 

{Pg - z)-^ = {Qa{z) + QA{xm + R{z) + Ra{z))-' = 0{l/h) 

Since e±^"'(^'''^) = 0{h-^l^+^)L2^L-^ for some M, it follows that 

{P-iW- z)-^ = C»(/i-^)i2^i2 z e [-5/2, 5/2] - i/i[0, (Ao - 3e)/2], , 
(P-iW^-z)-^ = 0(1/Imz)i2_,i2, Im2>0, 

where the second is immediate from non-negativity of 1^ as an operator. 

We now use a semiclassical maximum principle [7, Lemma 4. 7], [36, Lemma 2] to obtain 
the bound for (P — iW — z)~^ in (1.18) (after adjusting 5 and e). 

Remark 7.2. Strictly speaking we proved (1-18) for z e [—5/2, 5/2] — ih[Q, Ao/2 — eo], for 
any eo, provided that h is small enough. 

8. The CAP reduction of scattering problems: Proof of Theorem 3 

In this section we will prove a generalization of Theorem 3 which applies to a variety of 
scattering problems. Our approach of reduction to estimates for the Hamiltonian complex 

absorbing potential (CAP) is based on the work Datchev-Vasy [13] (see also [20, §4.1] and 
[15]) but as the argument is simple and elegant we reproduce it in our slightly modified 
setting. 

Let (y, g) be a complete Riemannian manifold and let 

(8.1) Pg = -h^Ag + V, V e C~(r; M). 

We make general assumption on {Y, g) which will allow asymptotically Euclidean and 
asymptotically hyperbolic infinities. 

Wc assume that Y is the interior of a compact manifold Y with a boundary, dY ^ 0. 
We choose a defining function of dY: 

(8.2) peC~(F;[0,oo)), {p = 0} = dY, dplay ^ 0. 
Let Pg — \^\l + V{x) be the principal symbol of Pg and let 

be the Hamiltonian flow (geodesic flow lifted to T*Y when V = 0. The flrst assumption 
on (Y, g) we make is a non-trapping (convexity) assumption near inflnity formulated using 
p with properties (8.2): 

(8.3) p(a;(t))e(0,eo), ^^p{x{t)) ^ ^ ^p{x{t)) < 0. 



(8.6) 
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The trapped set at energy E e [—S,S] is defined as 

(x,^) e Ke Pg{x,i) — E and exp(Mi7p^)(x, ^) is compact in T*Y. 
We assume that the trapped set at energies \E\ < 5, (see (1.15)), 

(8.4) is normally hyperbolic in the sense of (1.17) . 

We now make analytic assumptions on P. For that wc first assume that Pg can be 
modified inside a compact part of Y, to obtain an operator 

Poo = -h'Ag + v, ve c°°(r), vrp<.o= y\p<eo, 

with the foUowing properties: for some sq > and Cq > 0, 

(8.5) ||p^°(Poo -E- iO)-'p'°\\L2^Y)^L^Y) < ^, \E\ < 5, 

and 

u^{P^-E-iO)-'f, feC^iv) =^ 

WFh{u) \ WF,(/) C exp([0, oo)HnepJ (WF,(/) np^\E)) , 

where p^o = + V. 

We note that these assumptions do not require that the resolvent of Pqo has a meromorhic 
continuation from Im^; > to the lower half-plane. A stronger conclusion will be possible 
when we make that assumption: more precisely, for x ^ ^^{Y), we assume that the 
resolvent (Pqo — z)~^ continues from Imz > analytically to [—(5,(5] — i/i[0, Cq], for some 
Cq > 0, and that for some A^, the following resolvent estimate holds: 

(8.7) xiPoo - z)-\ = C»i,2^i2(/i-^), z e [-(5, S\ - %h% Co] . 

When Poo is chosen to be selfadjoint, interpolation [7, Lemma 4. 7], [36, Lemma 2] shows 
that (8.7) improves to 

(8.8) x(Poo-^)-'x = C»L^-.LK/i-'+^^^'"^/Mog(l//i)), [-(5,(5]-z/i[0,Co]. 
We can now state a more general version of Theorem 3: 

Theorem 6. Suppose that the Riemannian manifold (Y, g) and the potential V satisfy the 
assumptions (8.3), (8.4), (8.5) and (8.6). In particular, the trapped set for the operator 
P — —h^Ag + V is normally hyperbolic. 

Then, for some constant Ci (and Sq in (8.5) we have 

(8.9) \\p^^{Pg -E- iO)-V^1|L2(y)^L^(y) < Ci^^^^, \E\ < S. 

If in addition (8.7) holds then, for any e > 0, x{Pg " z)~^X can be continued analytically 
to [— (5/2, 5/2] — i/i[0, min(Co, Ao/2 — e)], with Aq given by (1.19), and 

(8.10) x{Pg-z)-\ = OL^^L<h-''), [-5/2,^2] -^/i[0,min(Co,Ao/2-e)], 
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with the improved estimate (1.18) if (8.8) holds. 

Before the proof we present two classes of manifolds to which satisfy our assumptions. 
We say {Y, g) is asymptotically Euclidean if 

g = p-^dp^ + p^'^goip), near dY, 

where go{p) is a family of metrics on dY depending smoothly on p up to p = 0. We say 
(y, g) is evenly asymptotically hyperbolic if 

g = p-'^dp'^ + p~'^go{p), near dY, 

where p is as before but this time go{p) is a family of metrics on dY depending smoothly 
on (hence even) up to p = 0. 

In both cases the non-trapping assumption near infinity (8.3) is valid: see [13, Proof 
of Lemma 4.1] for the asymptotically hyperbolic case; the asymptotically Euclidean case 
follows from the same proof, with the fourth displayed equation of the proof replaced by 
[38, (4.3)]. 

For asymptotically Euclidean manifolds (8.5) and (8.6) follow from the results of [38]. 
The modification of V can be done in any way which produces a non-trapping classical 
flow: for instance we can choose V = V + Vint where T^nt is a smooth, large non-negative 
potential (a barrier) supported in {p > eo}. 

To obtain (8.7) more care is needed but, under additional assumptions one can use 
an adaptation of the method of complex scaling of Aguilar-Combcs, Balslev- Combes and 
Simon - see [41] for the case of manifolds and for references. The simplest example for which 
this is valid was considered in Theorem 1. For even asymptotically hyperbolic manifolds 
the properties (8.5), (8.6), and (8.7) all follow from the recent work of Vasy [39]. 

As long we are not interested in analytic continuation properties, the weaker assumptions 
(8.5) and (8.6) may hold in the generality considered by Cardoso- Vodev [8]. 

Proof of Theorem 6. To show how Theorem 6 follows from Theorem 2 we use the parametrix 

construction of [13, §3]. For that we first have to relate the situation in this section to the 
set-up in Theorem 2. It will be convenient to rescale p so that in (8.3) we can take eo = 4. 

Let X be any compact manifold without boundary such that F C X is a smooth embed- 
ding: for example, we may take X to be the double of Y . We then extend p to p G L°°(X) 
to be identically on X \ y. Let P e ^'^{X) be any selfadjoint semiclassical differential 
operator satisfying 

P\p>l — Pg\p>l-, P = —h Ag^ + Vx, 
where gx is a Riemannian metric on X and Vx G C°°{X;M.). 
We then take W eC°°{X; [0, oo)) such that 



^ ^ [ 1 for p{x) < |. 
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Let V e C°°{Y) be a potential for which (8.5) and (8.6) hold. We notice that one possibility 
to obtain the required properties for Poo is to take a complex potential V — V — iW^ where, 
W^ooeC-(y;[0,oo)) 



for p(x) < 4, 

1 for p{x) > 5, 



see Fig. 3. 



Using the convexity property (8.3) it is easy to check that this operator satisfies (8.5) and 
(8.6). Then for Imz > 0, \ Rez\ < S, define the following holomorphic famihes of operators 



Rx{z) = {P-iW- z)-\ R^{z) = {P^^ 

Due to the compactness of X, the family of operators Rx{z): L'^{X) —f L'^{X) is mero- 
morphic for z E C. The resolvent Rx{z) is estimated in Theorem 2. For the moment we 
only assume that Roo{z) ■ L'^{Y) — >■ L'^{Y) is holomorphic for Im^; > and satisfies (8.5), 
(8.6). 

Now take a cutoff function xx e C°°(M, [0, 1]) with 

suppxx C (2, oo), supp(l - Xx) C (-00, 3). 
We put and Xoo = 1 - Xx- 



Xx{p+1) 



supp 1 1 




XX (p) \Xoo{p) 



Xoo(p- 1) 



STipp 



Figure 3. Schematic representation of the cut-offs used in the proof of The- 
orem 6 as functions of p{x). The spatial infinity is represented by p{x) — 
and X\Y corresponds to p{x) < 0. 

Our first Ansatz for the inverse of {Pg — z) is the operator 

F{z) = xxipi*) + l)Rxiz)xxipi*)) + Xoo(/5(«) - l)Rooiz)Xooipi*))- 

Note that F{z): L'^{Y) — )■ L'^{Y) for Imz > since all the cut-off functions are supported 
away from X \ Y. Also, the support properties of W, Woo and xx show that 

{Pg - z)xx{p{*) + 1) = Xxipi*) + l)iP-tW-z) + Ixxipi*) + 1), h'\], 

{Pg - z)Xoo{p{*) - 1) = XM*) - l){Poc -Z) + ixM*) - l),/i'A,]. 

Hence 

{Pg - z)F{z) = I + Ax{z)+ ^oo(^), 
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where 

Ax{z) = [xx{p{*) + l).h^^,]Rx{z)xx{p{*)), 

^ooW = [Xoo(p(«) - l),h^^,;^Roo{z)Xoo{p{*)). 

Note that Ax{zY = A^{zY = 0, due to the support properties 

supprf(xx(p(») + 1)) nsuppxx(p(»)) = 0, 
suppci(xoo(p(«) - 1)) nsuppxoo(p(»)) = 0- 
Moreover, thanks to assumptions (8.3) and (8.6) (see [13, Lemma 3.1]), 

(8.12) \\A^{z)Ax{z)\\l'^(y)-^l'^{y) = 0(/i°°), uniformly for Imz > 0, | Rez| < 5. 
Consequently 

(8.13) (P, - z)F{z){{I - Ax{z) - A^{z) + Ax{z)A^{z)) =1- E{z) , 

(8.14) where E{z) = A^{z)Ax{z) - A^{z)Ax{z)A^{z) . 

Using (8.12) we see that E{z) = 0{h°^)L2(^Y)^L'2{Y), uniformly for Imz > 0, |Re2;| < S. 
This allows to write an explicit expression for (P^ 



-1. 



9 



{Pg - zY' = F{z){I - Ax{z) - A^{z) + Ax{z)A^{z)) J] E^zf . 



n=0 



We now want to estimate ||p*°(Pg — z)~^ p^°\\l2^y)^l'^{y)- For this we expand the above 
identity using the expression of F{z) (some terms vanish due to the support properties 
(8.11)). Denoting ax = \\Rx{z)\\, a^o = \\p^°Roo{z)p^'^\\i we get the bound 

(8.15) \\p'°{Pg - z)-^p'° II <C{a^ + ax + ha^ax + h^l^ax) + 0{K^) . 

Finally, we use the bounds (8.5) for Coo, the bound (1.18) for ax (with Imz > 0), and 
obtain the desired estimate (8.9). 

When the assumption (8.7) holds, the construction shows that for x ^ C^(^) equal to 
1 on a sufficiently large set, 

oo 

x{Pg - z)-\ = X^Axij - Ax{z) - A^{z) + AxA^{z))xY.{E{z)xr. 

n=0 

continues analytically to the same region as both Rx{z) and xRoo{z)x- The same expansion 
as above allows to bound from above ||x(-P9 ~ -2)~^xll by the same expression as in (8.15), 
now using ax = ||xPx(2;)x||, Qoo = ||x^oo(^)x||- By using (1.18) for ax, resp. (8.7) for a^ 
(with now Imz taking negative values), we obtain (8.10). □ 

For completeness we conclude this section with the proof of Theorem 1. The conclusion 
is valid under more general assumptions of Theorem 6. 



56 S. NONNENMACHER AND M. ZWORSKI 

Proof of Theorem 1. In the notation of Theorem 6, (1.3) is equivalent to the estimate 
(8.16) M{P,)e-^'''^'WWiy)-,L^^y) < C^^^ e"^* + 0(/.~), 7 = ^(Ao - e) , 

vahd (with different constants) for any x ^ C^^X)- Let ip e C^(C) be an almost analytic 
extension of ■0, that is a function with the property that ^iid dz'i'iz) — 0{\ lm.z\'^) 

(sec for instance [43, Theorem 3.6]). We can construct ■0 so that supp^/; C [—5/2,5/2] — 
i[—5, 5]). We start with Stone's formula 

xV(P.)e-^*''«/'^X = ^ / V'(A)e-^*x((^.-A-iO)-^-(P,-A + iO)-^)xrfA. 

We now write R^{z) = {Pg — z)~^, for the resolvent in Im^; < (that is for the analytic 
continuation of {Pg — {z—iO))~^ from Im^ < 0) and R+{z) for the meromorphic continuation 
of the resolvent from Im 2; > to the lower half-plane. We then apply Green's formula to 
obtain, for < 7 < Aq, 

xHPg)e-''''^/'x = ^ / e-''^/\{R+{z) - R_{z))xi^{z)dz 

(8.17) 

+ - e-*^/\(i?+(z) - R4z))xdM^)dm{z) , 

J J-'yh<Imz<Q 

where dm{z) is the Lebesgue measure on C. Prom (1.18) (see Theorem 6) we get 
\\xR+{z)x\W^L^ < C/i-(^+^^)log(l//i), \\xR-{z)x\W^L^ < C/|Imz|, 

for — 7/i < Im^ < 0. Inserting these bounds in (8.17) gives (8.16) and that proves (a 
generalized version of) Theorem 1. □ 



9. Decay of correlations for contact Anosov flows: Proof of Theorem 5 

Most of this section is devoted to the proof of Theorem 4. This proof will be obtained by 
adapting the proof of Theorem 2, after reviewing the geometric point of view of Tsujii [37] 
and Faure-Sjostrand [21] (see also [12]). At the end of the section we deduce Theorem 5 
on the decay of correlations. 

9.1. Geometric structure. Let X be a smooth compact manifold of dimension d — 2k—l, 
k > 2. We assume that X is equipped with a contact 1-form a, that is, a form such that 
{da)^''~^ A a is non-degenerate. The Reeb vector field, S, is defined as the unique vector 
field satisfying 

e ker a^, axC^x) = 1, x e X . 

We assume that 



(9.1) 



'ft = exptS defines an Anosov flow on X. 
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That means that at each point x & X, the tangent space has a 7t-invariant decomposition 
into neutral (one dimensional), stable and unstable subspaces (each {k — 1) -dimensional): 

(9.2) ^Eoix)®E,{x)®Eu{x), Eo{x) ^RE^. 

We note that Eu{x) © Es{x) span the kernel of ker ax- 

The dual decomposition is obtained by taking Eq{x) to be the annihilator of Es{x) © 
Eu[x), E*{x) the annihillator of Eu{x) ®Eq[x), and similarly for E*{x). That makes E*[x) 
dual to Eu{x), E*{x) dual to Es{x), and Eq{x) dual to Eq{x). The fiber of the cotangent 
bundle then decomposes as 

(9.3) T;X = E*{x) © E;{x) © E:{x). 

The distributions E*{x) and El{x) have only Holder regularity, but Eq{x) and E*{x)®El{x) 
are smooth, and Eq{x) = M^a^ C T*X. 

The approach of [21] highlights the analogy between this dynamical setting and the 
scattering theory for the Schrodinger equation. The role of the Schrodinger operator is 
played by the (rescaled) generator of the flow — exp tE: 

(9.4) 7> = e''^/^'u, M e C°°(X), P = -ihE . 

The principal symbol of P simply reads p{x,^) = ^(S^.). 
The flow 7t can be lifted to a Hamiltonian flow on T*X: 

(fit : (x,0 ^ {jt{x)*djt{x)~^0, 
which is generated by p{x, ^): (fit — exjptHp. 

For each energy £■ e R, the energy shell p~^{E) is a union of affine subspaces: 

p-\E) = IJ {(x, : «.(0 = E}=\J {Ea, + El{x) + E^ix)) . 

We note that each of these energy shells has infinite volume; as opposed to the scattering 
theory setting, infinity occurs here in the momentum direction (the fibers), while the spatial 
direction is compact. 

The Anosov assumption implies that for t > 0, 

(9.5) \Mx,0\<Ce-''\^\, ^eE;{x), |(^_,(x, 01 < Ce-^ICI, ^ ^ Ki^), 

where | • | = | • denotes a norm on T*X, and we consider ifit{x,^) G T*^^^^^ ^^-^X . Since 
the action of fit inside each fiber T*X is linear, we see that the only trapped point in T*X 
must be on the line Eq{x). More precisely, the trapped set at energy E' G M is given by 

Ke^U {E*,{x) np-\E)) =\jEax, 

xex x&x 

that is Ke is the image of the section Ea in T*X. 
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Stacking together energies E e {1 — 6,1 + S), < 6 < 1, we obtain the trapped set 

= K= y Ke = {Ea^, X e X, \E -1\ < 6} C T*X . 

\E-1\<S 

This trapped set is normally hyperbolic in the sense of (1.17). 

Indeed, we first note that is a symplectic submanifold of T*X of dimension d+1 — 2k. 
Indeed, using (x, E), x E X, as coordinates on K^, {x, E) i->- Ea^, we have 

u\ks^ d{Ea) ^ dE Aa + Eda. 

This form is nondegenerate ior E ^ since a is a contact form. 

The tangent space to is given by the image of the differential of 

X xR3 {x,E) ^ Ea^ = (x, ^ = E/3{x)) , 
where we see (3{x) as the vector in M*^ representing a^- Hence, 

TEa^K' = E{daUT,X, .) + Ra, = {{v, E dp{x)v + s^{x)) ■.veT,X,se R} 
^ ' ^ CT^X® T*X = TEa. {T*X) . 

Here dP{x) can be interpreted as the Jacobian matrix d/3/dx on R'^"-'^. 

For each x e X, the symplectic orthogonal to Teu^^K^, denoted (TEa^K^)-^, can be 
obtained by lifting the vectors in ker as follows: 

V e ker a, ^ L^{v) = {v, EXd/3{x))v) e T,X © T^X = TEaAT*X) , 

where \d/3{x)) denotes the transpose of d(3{x). This subspace (TEa^K^)-^ is symplectic and 
transverse to K^: 

TpK^ ® {TpK^)^ = Tp{T*X), Vp = Ea^ e . 
Since kera^; = Eu{x) © Es{x), we can naturally split the orthogonal subspace into 

{TpKY = E; ® E; , E; = LUEu{x)), E; = LUEs{x)), p = Ea,eK'. 

The distributions E^^^ arc transverse to each other and fiow-invariant. E^^^ is a particular 
subspace of the global unstable subspace Eu{x) © -E*(a;) C TEax{T*X), and similarly for 
^Ea^ ■ Hence, in the present setting, the subspaces exactly correspond to the subspaces 
described in Lemma 4.1. 

9.2. Microlocally weighted spaces and the definition of resonances. Following [12] 
we now review the construction [21] of Hilbert spaces on which P — z (with P given in 
(9.4)) is a Predholm operator for Im^; > —^h, for some arbitrary ^ > 0. 

The key to the definition of these Hilbert spaces is a construction of a weight function 
which we quote from [21, Lemma 1.2] and [12, Lemma 3.1]. We use the notation E^ — 
U,^xE:{x)gT*X. 
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Lemma 9.1. Let Uo, U'q he conic neighbourhoods of E^, with Uq Uq and C/on(£'*U£'*) = 
0. There exist real-valued functions m e S^{T*X), f^ e S^{T*X) such that 

(1) m is positively homogeneous of degree for \$\ > 1/2, equal to —1,0,1 near the 
intersection of {\^\ > 1/2} with E*,Eq,E*, respectively, and 

(9.7) Hpm < near {Uq \ Uo) {\^\ > 1/2}, Hpm < on {\^\ > 1/2}; 

(2) (0~"^/o > c > for some constant c; 

(3) the function G = m log /o satisfies 

(9.8) HpG < -2 on > 1/2} \ Uo, HpG < on {\^\ > 1/2}. 

The function G also satisfies derivative bounds 

(9.9) G^O{\og{0), d:d^^H'^G^O{{0'\^\+^), |a| + |^| + ^>l, 
for any e > 0. 

As in [12, §3] we define 

(9.10) HtG{X)''^'e-''''"L\X,dx), 
where t > is a positive parameter. 

The domain of P acting on Hfc is defined as 

(9.11) VtG = {u e V'{X) : u,Pue Htc}. 

The action of P on Htc is equivalent to the action of the operator Pto on L^: 

P,a e*«™Pe-*^'" = exp(t adG»)P 

(9.12) fk 

= VTTad^.P + i?iv(a;,/iD), Rn e Ve > 0. 

k=0 

The vahdity of (9.12) follows from the fact that the operators e^^'^'" are pseudodifferen- 
tial operators [43, Theorem 8.6], hence the pseudodifferential calculus apphes directly [43, 
Theorem 9.5, Theorem 14.1]. This expansion and the arguments in [21, §3] give 

Proposition 9.2. For Pto defined by (9.12), we have 

i) the operator Pto — z : T>{PtG) — >■ L"^ is Fredholm of index zero for Im^; > —th. Here 
V{PtG) is the domain of Pto- 

a) for C > large enough, {Pto — z) is invertible on {Imz > Ch}. 

In [21] the above construction was performed, replacing the /i-quantization by the h — 1 
quantization. It lead to the construction of HtG,i{X) = e~*'^"'f^'^^L^(X) equal, as vec- 
tor space, to the above /i-dcpendent space HtQ{X). The norms of these two spaces are 
equivalent with one another, but in an /?.-dcpcndcnt way: 

(9-13) h^\\A\HtG,^{x)/Co < \\u\\HtG < Coh~^ ll^ll/ZtciW > 
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- see [12, §3]. As a consequence, if we call Pi — — iS, Theorem 4 translates into the fact 
that PtG,i = e^G-'{x,D)p^^-tG-'{x,D) -g fvedholm in the strip {ImA > -t}, admits finitely 
many eigenvalues in the strip {ImA > — A/2 + e}, and satisfies the resolvent estimate 

(9.14) ||(PtG,i-A)-i^2^i2 = 0(A^°), ImA>-A/2 + e, |ReA|>C. 

9.3. Reduction to Theorem 2. In order to prove Theorem 4, we proceed as in the proof 
of Theorem 6 in §8, by constructing two operators which microlocally agree with Ptc (up 
to negligible error terms) on different subsets of T*X. 

Let W e ^^(X) be as in Example 2 of §1.3. The trapped set defined in §1.3 agrees 
with the trapped set in §9.1 and, as shown in (9.5), it satisfies, the assumptions of normal 
hyperbohcity. Hence Theorem 2 apphes to P = P - iW. If A e ^^°™p(X) satisfies 



We now introduce an operator which has better global properties and agrees with Ptc 
near infinity. For that we proceed as in the proof of Theorem 6, and take Woo G ^'^°™p(X) 
such that 



(9.17) WFhiWoo) <^{G-\0)y, WFh{I -Woo) c{:k', WFh{W) nWFh{Woo) = 0. 



Properties of the operator Poo are listed in the following 

Lemma 9.3. Fix (3 > and let t be large enough so that Pta — z is a Fredholm opeartor 
for Imz > —^h. Then, there exists Nq and Hq such that, for < h < Hq, 



Proof. The first part follows from the now standard non-trapping estimates (see [35, §4]). 
In the setting of Anosov flows the details are presented in the proof of [12, Lemma 5.1] 
(only the escape function constructed in Lemma 4.5 above is needed). 





Then for any B with WF^J - P) C ZWFh{Woo), 

(9.18) (/ - P)Poo = (/ - B)PtG + 0{h'^)v'^o 
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The propagation result is a real principal type propagation result [43, Theorem 12.5] 
which holds when the imaginary part of the symbol is non-positive - see Lemma A.l below 
for a dynamical version. □ 




|e|->oo 

Figure 4. Schematic representation of pseudodifferential cut-offs used in the 
proof of Theorem 4. The horizontal axis corresponds to |^|, the cotangent 
variable. Infinity in |^| plays the role of p = in Fig. 3. The asymmetry 
is intentional, to stress that there is no need for an auxiliary manifold, as 
opposed to the proof of Theorem 6 illustrated in Fig. 3. 



Proof of Theorem 4- The proof is a repetitition of the proof of Theorem 3 with Rx re- 
placed by (P — z)'^ and Roo by (Poo — The spatial cut-off functions are replaced by 
pseudifferential operators: xx{p{x)) is replaced by Aq G ^^°™'^{X), satisfying 

WFhiAo) n {{x, : leis. > M} = 0, WF^(7 - ^o) n WF^iW^) = 0, 

where M is given in the definition of W, see (1.13). The function Xx{p{x) + 1) is replaced 
by Ai e ^'^°"'P(X), where 

WF;,(/ - A,) n WF,(ylo) = 0, WF;,(yli) n {{x, : \^\g^ >M} = 0, 

xMx)) is replaced by Pq = I - Ao E ^'°(X), and finaUy xMx) - 1) by Pi e *°(X) , 
where 

WF^(w^oo) n WF^(Pi) = 0, WF;,(/ - Pi) n wf^Pq) = 0. 

We also require that 

WF^(Ai),WF^/-Pi) c (6-^(0))°. 
The parametrix is now obtained by putting 

F{z) = A,{P-tW- z)-'Ao + B,{PtG - tWoo - z)-'Bo. 

Using (9.16), (9.18) and Lemma 9.3 we obtain the theorem by proceeding as in the proof 
of Theorem 3 in 58. □ 
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Proof of Theorem 5. We will use the nonsemiclassical operator Pi — — iS. It is selfadjoint 
on LF'{X) - see [21, Appendix A] - hence, by Stone's formula, we get for any f,g& C°°{X) 

f ^*_Jgdx={e-''''^f,g) 
Jx 

= I e-^'* {m - A - iO)- V, -g) - {{Pi - A + iO)- V, -g)) dX 

= ^Et / e-'^\^ + i)-''{{Pi-X±iO)-\Pi + i)'f,g)dX. 

Here the brackets (•,•) represent L'^{X) scalar products. 

For t > we can deform the contour in the integral corresponding to +iO (A approaching 
the real axis from below), where ||(Pi — A)~^|| < | Im A|~^, so that for /c > 1 the integral is 
bounded as 

-7^ [ e-^*(A + t)-'{{Pi - X)-\Pi + iff, -g)dX = 0(e-*^||/||^.||y|U.) . 

Thus, 

I Ytfgdx^^. I e-^'\X + i)-'{{Pi-X-iO)-\Pi + iff,-g) + Oae-'^), 

for any A, with the bounds depending on seminorms of / and g in C°°. We now use the 
nonsemiclassical weights G'^{x, D) constructed in §9.2 to conjugate Pi, and write 

/ ilfgdx^ 

Jx 

-i. / e-^'\X + i)-'{{PtG,i - A - iO)-^(P,G,i + i)^e*«"'(^'^)/, e-'^^^^'^^g) + 0/,,(e-*^). 

The nonsemiclassical analogue (9.14) of Theorem 4 shows that, by taking A; > A'^g + 1, we 
may deform the contour of integration down to Im A = — Ao/2 + e, collecting finitely many 
poles jij, to finally obtain the expansion (1.22). □ 

Appendix: Evolution for the CAP-modified Hamiltonian. 

In the appendix we show some properties of the CAP-modified Hamiltonian, that is the 
Hamiltonian modified by adding a complex absorbing potential. At first we work under 
the general assumptions (1.9). 

The semigroup exp(— it(P — iW)/h) : L'^{X) — )• L'^{X) is defined using the Hille-Yosida 
theorem: for h small P — iW — i is invertible as its symbol is elliptic in the semiclassical 
sense (see (1-11) and [43, Theorem 4.29]). Ellipticity assumption for large values of ^ also 
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shows that P — iW is a Predholm operator, and the comment about invertibihty shows 
that it has index 0. The estimate 

||(P-iiy-^)ii||||ii|| > -Im{(P -iW - z)u,u) > Im^ll-ull^ u E HJ^^'iX), 

then shows invertibihty for Imz > 0, with the bound 

II {P-zW- z)-^ \\l2^l2 <v^, lmz>0. 

Imz 

Since the domain of P — iW is given by H"^(X) which is dense in L^, the hypotheses of 
the Hille-Yosida theorem are satisfied, and 

||e-^*(^-^^)/'^|U.^^. < 1, t>0, 

^ ■ ^ ^-it{P-iW)/h^-is{P-iW)/h ^ ^-i{t+s){P-iW)/h^ t,s>0. 

Alternatively we can show the existence of the semigroup exp{—it{P — iW)/h) using energy 
estimates, just as is done in the proof of [43, Theorem 10.3]. We get that for any T > 0, 

(A.2) e-itiP-iW)/h g c{[0,nC{HfXX),Hf^{X)))nC\[0,T]-Cim,K'n)- 

Our final estimates will all be given for only and that is sufficient for our purposes. 

The first result we state concerns propagation of semiclassical wave front sets. We recall 
the notation (pt — exp{tHp) for the Hamiltonian flow generated by p{x, ^). 

Lemma A.l. Suppose that A G Then for any T independent of h there exists 

a smooth family of operators 

(A.3) [0, T]3t^ Q{t) e *^°'"P(X), WF;,(/ - Q{t)) n ipt{WFh{A)) = 0, 

such that 

(A.4) (/ - Q{t)) e-**(^-^^)/M = 0{h^)L2^L2. 

In addition ifWFh{A) C w~^{[e, oo)), e > 0, then for any fixed t > 0, 

(A.5) ^-it{P-iW)/h^ ^ C»(/i~)i2_^^2, A e-^*(^-^^)/'^ = 0{h'^)L2-,L2. 

Proof. We first construct Q{t) using a semiclassical adaptation of a standard microlocal 
procedure - see [27, §23.1]. For that, let Q{0) e *^°'^p(X) be an operator satisfying 
WF h{I — Q{0)) n WF h{A) — 0, and with the principal symbol, qo{0), independent of h. 

Using the fact that the fiow Lpt is defined for all t we put qo{t) =^ Lp*^^qo{0). In terms of the 
Poisson bracket on the extended phase space T*(]Rt x X) 3 (t,x,r,^), this means that the 
function qo{t) satisfies the identity {r +p, ?o(^)} = 0. Consequently, at the quantum level 
we have 

[hDt + P, Opngo(t)] = hR^it), R,{t) e *--p(X), 
Opngo(O)) - g(0) = hE„ E, e *--p(x). 
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and the principal symbols of Ri, Ei, ri,ei e C'^{T*X), are independent of h. If pi = 
(t((P — Op)^(p))//i, we then solve (in the unknown qi{t)) the equation 

{T+p,gi(i)} = ri - {pi,?o(*)}, 9i(0) = ei. 
By iteration of this procedure we obtain G C°°(T*X) such that 

AT-l 

[/lA + h'Ovliqm] = h'^RNit), RN{t) G 

£=0 

7V-1 

E h^Op^Qom - Q{0) = /i^£;jv, Em g ^^"'^^(X). 
By a standard Borel resummation we may construct Q{t) G ^'^°™p(X) such that Q{t) ~ 

E,>o/^^Opn*(t)). 

For any iV > we can iteratively construct a sequence of auxiliary operators Qj{t) — 
Qj{t)* G *"°'"P(X), < J < A^, satisfying 

WF,,(/ - g,+i(t)) n WFhiQM = WFh{I - QM n ^t{'^Yn{A)) 
(A.6) = WF;,(/ - Q{t)) n WYn{Qj{t)) = 0, 

[Q^ (t), /lA + P] G ([0, T]; /?,°°^'=°'^P(X)) . 

(These assumptions imply that ^pt{W¥h{A)) C WF(Qj(t)) C W¥{Qj+i{t)) C WF(Q(i)).) 

= 1. Our aim is to prove the following property: 



(A.7) Wj{t) = (J - Qj{t))v{t) = 0{y''')L2, for J = 0, . . . , iV, < t < T. 

Since A G ^P^^P, (A. 2) shows that this property holds for j — 0. Let us now prove that, if 
true at the level j, it then holds at the level j + 1. 

Noting that 

(A.8) = (/ - Qj+i)wj + C»(/i~)c-, 

we have 

{hDt + P- iW)w,+i = {I- Q,+i{t)){hDt + P- iW)wj - ^[l^,Q,+i]u^,- + 0{K^)l2 
Dividing hy h/i, taking the inner product with Wj+i, taking real parts and integrating gives 



hj+lWllL2 



ki+i(0)||i2 +2 j {Wwj+i{s),w^+i{s))ds^ 

(A.9) ^ 

^ Re([W^,Q,+i(s)]w,(s),w,+i(s))ds + 0(/i°°), 



h 
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Now, 

(/ - Qj+i{s))[W, (/ - Q,+i(s))] = thBj+,{s) + h^C,+,{s), 

Hence, using (A. 8) and the induction hypothesis (A. 7), the right hand side of (A. 9) becomes 

2h I Re{Cj+i(s)wj(s),Wj(s))ds + 0(h°°) ^0(h^+^). 
Jo 

Returning to (A. 9) and using the non-negativity of W, we see that 
Since 

«;,+i(0) = (/ - Q,+i)Au = 0{h^)L2, 
we have estabhshed (A. 7) with j replaced by j + 1. 
The estimate (A. 4) then follows from 

(/ - Q{t))v{t) = (/ - Q{t))wj{t) + OL2{h'^). 
the estimate (A. 7) at the level j = N, and the fact that N could be taken arbitrary large. 
To see (A.5) we note that if A e *^°°^p(X) then 

WFh(A) C w-\[e,oo) =^ ipt(WFh(A)) C w-\[e/2,oo) iorO<t<5. 
Hence, by (A.4), 

WFh{v{t)) C w-\[e/2, oo), v{t) = e~''^''-'^^/''Au, \\u\\l^ = 1, 0<t<5. 
This means that we can modify W into Wi, so that 

a{Wi){x, > (OVC, W, > Co, for < /i< ho, 

while we have 

= (hDt + P- iW)v{t) = {hDt + P- iWi)v{t) + O{h°°)coo uniformly ior < t < 5 . 
Taking the imaginary part of the inner product of the above expression with v{t) gives 

and hence 

\\v{t)\\l2 = Oih"^) uniformly for 6/2<t<6. 

This proves the first part of (A.5). The second part follows by taking a conjugate: 
j^^-tt{p-tW)/h _ ^^-tt{-p-tW)/hj^*Y ^ ^j^^ the arguments remain valid for P replaced 
by -P. □ 

The next lemma is needed in §7 and follows immediately from Lemma A.l: 
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Proposition A.2. Suppose that A e ^f^°°^P(X) satisfies 

(A.IO) WF,,{A) c p~\{-S,S)) nw-\[0,eo)), 

for some eo > and that T is independent of h. 

Then there exists B e ^'^°"ip(X) for which (A.IO) holds with B in place of A, and 

(A.ll) ^-iKP-iW)/hj^ ^ Be-'*^^-'^^>^A + 0{h'^)L2^L2, 0<t<T. 

Proof. Using again the operator Q{t) constructed in the proof of Lemma A.l, we take 
a compact set L containing WF h{Q{t)) for all < i < T. By taking WF/,(g(0)) C 
p~^{{—S,5)) (which is possible due the assumptions on A) we see that we can assume 
L C p~^{{-5, 5)). We can now choose B e ^''=°™p(X) such that 

WFft(7 - S) n L n w-\[0, eo/3]) = 0, WF^(S) C p-\{-S, S)) n w-\[0, eo/2). 

This implies that (/ - B)Q{t) = C(t), where WF,j(C(t)) C w-^([eo/3, oo)), and hence, by 
(A.4) and (A.5), 

(/ - 5)e-'*(^-'^)/'^/l = {C{t) + {I-B){I- g(t))e-**(^-*'^)/M = C(/i°°)L2^i2, 
proving (A.ll). □ 

Finally we present a modification of [32, Lemma A.l]. The modification lies in slightly 
different assumptions on P and W, and the proof also corrects a mistake in the proof given 
in [32]. From now on we work under the extra assumption (1.10) on the CAP. Wc remark 
that in [32] we only needed Lemma A.l and hence the assumption (1.10) was not required. 

Proposition A. 3. Suppose that X is a compact manifold, P is a self-adjoint operator, 
P G ^'"(X), W G ^*^(X), W > 0, and that (1.9) and (1.10) hold. Then for any t 
independent of h, for A G \[''=°™p(X) satisfying (A.IO), we may write 



where 



(A.12) 



^itP/h^-itiP-iW)/hj^ = VA{t) + 0{h'^)L2^L^, 

VA{t)e^r^{x), WF,{VA{t))c f]{^_,{w-\o)))nWF,{A), 

0<s<t 

a{VA{t)) = exp (^-^ iflWds^ a{A) . 



The class of operators was introduced in ^3.2. 

The proof is based on the following lemma inspired by the pseudodifferential approach 
to constructing parametrices for parabolic equations presented in 
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Lemma A.4. Suppose that 1 1->- p{t, z, h), p{t, •, h) e C^(R^";R), is a family of functions 
satisfying 

d!^d:p{t,Z,h)^Ok,a{i), P>-Ch, 0<h<ho, 

\d:p{t,z,h)\=0.{p'-'), 0<5<l 

Then, for < s < t there exists E{t, s) G \E'5(]R"') such that 

{hdt+p'"{t,x,hD^,h))E{t,s) ^0, t>s>0, E{s,s) ^ I. 

Moreover, E{t,s) = e^{t,s,x,hDx.ih) where e{t,s) e ^'^(R^") has an explicit expansion 
given in (A. 26) below. 

Proof Replacing phyp+{C + l)h, gives p > /i and p{t, ; h) e {C + l)h + C~(R2"). The 
multiplicative factor e^'^+^^^*~'^^ in the evolution equation is irrelevant to our estimates. 

For any A" > we try to approximate the symbol e(i, s, ^, h) by an expansion of the 
form 

N 

(A. 14) /iv(t, s, z, h) = ^ yej{t, s, z, h) . 

3=0 

The symbol of the operator hdtffi f^ can be expanded using the standard notation 
o = {a^b)"^ and the product formula (see for instance [43, Theorem 4.12]): 

hdtfN{t,s) + \p{t)#fN{t,s)] 

= (hdtej{t,s) + XI ^ [k^huj{D„,Dy,)Y p{t,z)ej{t,s,w%=^ + h^'^N, 

j=0 k=0 

J2 [{hdt + p{t))ej{t, s) + Y, ilMD,, D^))''' p{t, z)ee{t, s, «;) 



N-l 



+ h^rnit, s, z), rN{t, s, z) ^ rN,j{t, s, z). 

3=0 

The remainders satisfy the following bounds (see for instance [35, (3.12)]): 

sup |(9"rwj(t,s,z)| < 

z 

^^■^^^ Co.,N,j V sup sup {h^d,,^f{a{D,,D^)Y-^d'i^p{z)dl-e^{w) 

ai-\-a2=o: \^\— ^^^^ 

The standard strategy is now to iteratively construct the symbols ej so that each term in 
the above expansion vanishes. The term j = simply reads {hdt + p)eo = 0. From the 
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initial condition eo{s,s) = 1, it is solved by 

(A. 16) eo{t, s, z, h) = exp J^p{s', z, h)ds'^ . 

For j > 1, the symbol Cj is obtained iteratively by solving 

ej{t,s,z) = ^ j\o{t,s',z)qj{s',s,z)ds', ej{t,s,.) G ^(M^"), 

(A.17) 

qj{t,s,z) — — {'^tuj{D,,D^)Y~'p{t,z)ee{t,s,w)U=^ G C^i^T) ■ 

This construction formally leads to an approximate solution: 

(A.18) hdtfN{t, s, z) + [p{t, •)i^fNit, s, •)] (z) = h^r-Nit, s, z). 

To make the approximation effective, we now need to check that the sum (A. 14) is indeed 
an expansion in power of h. We thus need to estimate the e^'s and thereby the remainders 
rATj's. 

We will prove the following estimate by induction: 
(A.19) \d':e,{t,s,z)\ < Ca,jh-''^-'\"\ (^1 + ^)d8')''^'"') eo{t,s,z). 

For that we first note that, as p > h, and \d"p\ < CaP^~^, we have 
(A.20) < C„/i-^p. 

Consequently, for j = we have 

k 

\d:eo{t,s,z)\ < Yl ll{T.fs\9'''Pi^',^)\ds')eo{t,s,z) 
(A-21) <c„ J2 fl{h-'Ulp(''^^)ds')eo(t,s,z) 



< C>-^l"l (^1 + z)ds') eo(t, s, z), 



Here we used the fact that A; < |q;| and that 

A'' < c«(l + ^1-1), A = i z)ds' > 0. 

This gives (A.19) for j = 0. 

To proceed with the induction we put 

aj^a{t, s, z) = d"ej{t, s, z)/eo{t, s, z), bj,a{t, s, z) =^ d'^qj{t, s, z)/eo{t, s, z), 
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noting that, for some coefficients, c,, 

(A.22) ''^+^^=" 

aj,a{t,s,z) = - cp^AJ ao,f)At^s\z)hj,p^{s',s)ds\ 

wfiere the last lequahty follows from eo(i, s', 2;)eo(s', s, z) — eo{t, s, z), s < s' < t. 
Our aim is to show 

(A.23) \b,,^{t, s, z)\ < C^,,h-''^-'\Mt, z) (^1 + , 

and 

(A.24) \a,At, s, z)\ < C'^,h-''^-'\-\ (^1 + (i />(.', , 

assuming the statements are true for j replaced by smaller values. 

We note that the case of j = has been shown in (A.21), and since bo,a = 0. 

The first estimate (A.23) follows immediately from the inductive hypothesis on ae^a, 
< ^ < J — 1 and the estimates on p in (A. 20). The second estimate (A.24) follows from 
(A.21), (A.23) and the obvious fact that J^^' p{s')ds' < f^p{s')ds', s < Si < S2 < t. 

We note that (A. 19) and the definition of Cq given in (A. 16) imply that 

d^ejit, s, z) = C(/i-^l"l-2'5j) , j > . 

so from (A. 14) we see that the symbol f^it, s) e ^'^(R^"). 
The bounds (A. 15) then show that the remainders satisfy 

Going back to (A. 18) we get the expression 

(A.25) E{t, s) = f^{t, s, X, hD,) + h""-' ^ E{t, s')r]^(s', s, x, hD,). 

(We note that, since p'^{t,x,hDx) > —Ch by the sharp Garding inequality [43, Theorem 
4.32], and since p^ is bounded on L^, the operator E{t,s) exists and is bounded on L^, 
uniformly in /i.) Since operators in ^-re uniformly bounded on LP' [43, Theorem 4.23], it 
follows that 

Bit, s) = f^{t, s, X, hD,) + 0(/i(^-2^)^)i2^^2. 

To show that E{t,s) - e^{s,t,x,hD^) G ^^°"'P(M"), we use (A.25) and Beals's lemma in 
the form given in [35, Lemma 3.5, h = 1]: £j are linear functions on M^", £J = £J{x, hD), 
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then 

adiw ■ ■ ■ adj^w E{s, t) — 

ade^ ■ ■ ■ a.dq f^{s, t, x, hD^) + h^~^ J ad^™ • • • ad^^ {E{s, s')r'%{s' , s, x, hD^)) ds' 

if is large enough. Here we used the fact that fiq, £ Ss and that ad^y • • • ad^^- E{s, t) — 
0{1)l2^x.2, which follows from considering the evolutions equation for the operator on the 
left hand side. 

In conclusion we have shown that E{t,s) — e'^{t, s,x, hD^), where e e ^'^(M") admits 
the expansion 

(A.26) ^(^' ^' ^' ^) ^ E ^'^^•(^' ^' ^' ^^•(^' ^) ^ J > i> 

with Co given by (A. 16). □ 

Proof of Proposition A. 3. We first observe that Lemma A.l (applied both to propagators 
hrP-iW and for P) shows that for B e ^^°"^p(X) satisfying WFh{I - B)nWFh{A) = 0, 

We can choose B = B* . Since 
it follows that 

(A.27) B ^tPlh^-it{P-iW)ihj^ ^ yB(^) ^ 

where 

(A.28) /latV^^(t) = -l^B(t)V^''(t), W^ij(t) = Se'*-P/'^iye-'*-^/'^S. 

We note that Wb^I) e *'=°'^p(A:), WFh{WB{t)) C WF;,(S), and that Wsit) > 0. Hence 
V^{t) = C(l)L2_j.i2 and (A.27) follows from Duhamel's formula. 

By decomposing >1 as a sum of operators, we can assume that WF h{A) is supported in 
a neighbourhood of a fiber of a point in X. Hence, by choosing B with a sufficiently small 
wave front set, we only need to prove that V^{t) G ^5 for X = M"; that follows from 
Lemma A. 4, since the symbol of Wsit) satisfies the assumptions (A. 13). The second and 
third properties in (A. 12) follows from (A. 5) and (A.26). □ 
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